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ABSTRACT 

We investigate signs of Active Galactic Nucleus (AGN) in the luminous infrared galaxy NGC 3256 
at both infrared and X-ray wavelengths. NGC 3256 has double, the Northern and Southern, nuclei 
(hereafter, N and S nuclei, respectively). We show that the Spitzer IRAC colors extracted at the S 
nucleus are AGN-like, and the Spitzer IRS spectrum is bluer at < 6 ftm than at the N nucleus. We built 
for the S nucleus an AGN-starburst composite model with a heavily absorbed AGN to successfully 
reproduce not only the IRAC and IRS specrophotometries at ~ 3" but also the very deep silicate 
9.7 /im absorption observed at 0'.'36 scale by Diaz-Santos et al. We found a 2.2 pm compact source at 
the S nucleus in a H5TNICMOS image and identified its unresolved core (at 0'.'26 resolution) with the 
compact core in previous mid-infrared observations at comparable resolution. The flux of the 2.2 pm 
core is consistent with our AGN spectral energy distribution model. We also analyzed a deeper than 
ever Chandra X-ray spectrum of the unresolved (at 0'.'5 resolution) source at the S nucleus. We found 
that a dual-component power-law model (for primary and scattered ones) fits an apparently very hard 
spectrum with a moderately large absorption on the primary component. Together with a limit on 
equivalent width of a fluorescent Fe-K emission line at 6.4 keV, the X-ray spectrum is consistent with 
a typical Gompton-thin Seyfert 2. We therefore suggest that the S nucleus hosts a heavily absorbed 
low-luminosity AGN. 

Subject headings: infrared: galaxies — X-rays: galaxies — galaxies: active — galaxies: individual 
(NGC 3256) — galaxies: nuclei 

1. INTRODUCTION 

Occurrence of Active Galactic Nuclei (AGNs) in Lu¬ 
minous Infrared Galaxies (LIRGs) has been extensively 
discussed since there are many lines of evidence for direct 
and causal relationship between them. LIRGs predomi¬ 
nantly emit at far infrared (FIR), and the infrared (IR) 
luminosities at 8-1000 pm, Air, amount to 11 < log Air 
(Lq) < 12 by definition (Sanders & Mirabel 1996). Al¬ 
though bulk of the energy in LIRGs is due to intense star 
formation activity (Sanders & Mirabel 1996), AGNs are 
often found among LIRGs. Also, a significant fraction 
of LIRGs are tidally interacting objects. The fraction of 
AGN-dominated sources is known to be higher in LIRGs 
with larger IR luminosity (e.g., Veilleux et al. 1995; Kim 
et al. 1995; Nardini et al. 2008, 2010; Petrie et al. 2011; 

Iwasawa et al. 2011; Alonso-Herrero et al. 2012) and in 
later interaction stages (Petrie et al. 2011; Iwasawa et al. 

2011 ). 

Both mid-infrared (MIR) and X-ray studies are very 
effective to reveal nuclear activities in LIRGs. Although 
LIRGs, in particular their nuclei, are often heavily ab¬ 
sorbed at optical and near-infrared (NIR) wavelengths, 
both MIR and hard X-ray photons can penetrate through 
such interstellar medium. We can study AGN indica¬ 
tions at MIR because AGNs typically show character¬ 
istic power-law-like spectral energy distribution (SED), 
which often shows an excess at MIR over the SED from 
warm and cold dusts (e.g., Granato et al. 1997; Alonso- 
Herrero et al. 2001; Nenkova et al. 2002; Armus et al. 

2007; Donley et al. 2012). Various kinds of AGN di¬ 
agnostics have been developed at MIR, and empirical 
assessments of relative AGN contribution to the entire 


system have been made (Genzel et al. 1998; Sturm et al. 
2002; Armus et al. 2007; Spoon et al. 2007; Desai et al. 
2007; Petrie et al. 2011; Stierwalt et al. 2013). We can 
also study typical AGNs at X-ray because AGNs show 
luminous and characteristic power-law spectra, the hard 
part of which suffers from less absorption by circum-AGN 
material. Therefore, X-ray has been heavily utilized to 
study AGN in LIRGs (e.g., Iwasawa et al. 2011). 

NGC 3256 is the most luminous LIRG in the local uni¬ 
verse {z < 0.01). At a distance of 35 Mpe (for consistency 
with Sakamoto et al. 2014, which adopted Sanders et al. 
2003^), its IR luminosity (Air) is as large as 3.6 x 10^^ 
Lq (Sanders et al. 2003), or total IR luminosity at 3- 
1100 pm (TIR) is (2.7-3.3) x 10^^ Lq(E ngelbracht et al. 
2008 converted for our assumed distance). It is also 
among the most luminous X-ray sources without a con¬ 
firmed AGN in the local universe (Moran et al. 1999; 
Lira et al. 2002; Pereira-Santaella et al. 2011). It is a ma¬ 
jor merging system showing tidally distorted morphology 
at galaxy scale with tidal tails, disturbed outer spiral 
arms, and prominent dust lanes with complex morphol¬ 
ogy (e.g., Graham et al. 1984; Kotilainen et al. 1996; 
Zepf et al. 1999; Lipari et al. 2000; Alonso-Herrero et al. 
2002; English et al. 2003; Alonso-Herrero et al. 2006a,b; 
Sakamoto et al. 2014). Its double nuclei, often referred 
to as Northern and Southern nuclei (hereafter, N and S 
nuclei, respectively) are separated by only 5" (850 pc). 
They are clearly visible at both radio and MIR (e.g., Nor- 

^ Sanders et al. (2003) adopted distance with cz using the cosmic 
attractor model outlined in Appendix A of Mould et al. (2000), 
using Hq = 75 km s~^ Mpc“^ and adopting a flat cosmology in 
which = 0.3 and = 0.7. 
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ris & Forbes 1995; Kotilainen et al. 1996; Neff et al. 2003; 
Alonso-Herrero et al. 2006b; Lira et al. 2008; Sakamoto 
et al. 2014), but the S nucleus is hidden in dust at the 
optical. The merger is likely in the late stage just before 
coalescence of the nuclei (Armus et al. 2009; Stierwalt et 
al. 2013). 

The N nucleus is a core of a starburst galaxy. Its 
optical spectrum shows H ii region-like features. Ex¬ 
tended outflowing ionized gas with shocks is also de¬ 
tected and attributed to a superwind powered by the 
starburst (Moran et al. 1999; Lipari et al. 2000; Rich et 
al. 2012; Bellocchi et al. 2013). Stellar population ex¬ 
amined in K band indicates young starburst population 
(Doyon et al. 1994). Prominent Polycyclic Aromatic Hy¬ 
drocarbon (PAH) features as well as low-ionization fine 
structure lines have been detected at MIR (Siebenmor- 
gen et al. 2004; Martin-Hernandez et al. 2006; Lira et 
al. 2008; Pereira-Santaella et al. 2010; Diaz-Santos et al. 
2010a), also indicating star formation activity. The nu¬ 
cleus is spatially resolved with HST at 1.6 /xm at 0"14 
full-width half maximum (FWHM) resolution (Alonso- 
Herrero et al. 2002), the 8.1 m Gemini Telescope at 
8.74 /im at its diffraction-limit (O'.'SO FWHM) resolution 
(Alonso-Herrero et al. 2006b), and Chandra at X-ray at 
0'.'5 FWHM resolution (Lira et al. 2002). 

Presence of an AGN at the optically obscured S nu¬ 
cleus has been suspected for a long time (e.g., Kotilainen 
et al. 1996), but there is still no firm evidence for it. Dust 
lanes cover the S nucleus and make it invisible at optical 
and NIR wavelengths below K band. A compact (unre¬ 
solved at 0'.'5 FWHM) and moderately absorbed (neutral 
Hydrogen column density A^h — 5 x 10^^ cm“^) X-ray 
source was detected with Chandra at this nucleus (Lira 
et al. 2002). Although it looks much fainter (by at least 
two orders of magnitude) than expected for a classical 
Seyfert nucleus given its MIR luminosity, it is brighter 
than expected for typical starburst galaxies (Lira et al. 
2002; see also Awaki et al. 1991; Turner et al. 1997; 
Guainazzi et al. 2005; Fukazawa et al. 2011). Analysis 
of the MIR SED extracted within a kpc-scale aperture 
(including both nuclei and the host galaxy component) 
indicated that AGN contribution to the total MIR lumi¬ 
nosity is < 5% (Alonso-Herrero et al. 2012). Sakamoto 
et al. (2014) recently discovered a bipolar collimated jet¬ 
like molecular-gas outflow from the S nucleus, and argued 
that it is likely driven by an AGN. 

Due to its proximity, NGC 3256 is among the best tar¬ 
gets to examine the characteristics and roles of AGN, 
if any, within a LIRG system. For this purpose, it is 
essential to verify the presence of an AGN in the S nu¬ 
cleus using spatially resolved and/or sensitive MIR and 
X-ray observations (e.g., Lira et al. 2002; Diaz-Santos et 
al. 2010b, 2011; Pereira-Santaella et al. 2010). Recent 
spatially-resolved MIR spectroscopies revealed some in¬ 
teresting differences between the two nuclei. By using the 
slit-scanning spectroscopy data taken with the Spitzer 
IRS (Houck et al. 2004), Pereira-Santaella et al. (2010) 
noted that the S nucleus shows deeper silicate absorp¬ 
tion at 9.7 /xm, stronger H 2 lines, and a slightly smaller 
equivalent width (EW) of PAH 6.2 /xm than those at 
the N nucleus. More dramatic differences are found in 
ground-based spectra obtained with higher-spatial reso¬ 
lution. The silicate absorption is much deeper and the 


PAH features are undetectably weaker at 0736 (61 pc) 
scale at the S nucleus than in the IRS spectrum (Dfaz- 
Santos et al. 2010a). At X-ray, Lira et al. (2002) reported 
that the point-like source at the S nucleus is more heavily 
absorbed than other point sources and the diffuse com¬ 
ponent within the galaxy. In this work we utilize archival 
and published data of Spitzer, HST, and Chandra to an¬ 
alyze the nuclear spectrophotometry more comprehen¬ 
sively than ever to explore an AGN in the S nucleus. 

2. INFRARED DATA 
2.1. Spitzer IRAC 

We examined Spitzer IRAC (Eazio et al. 2004) archival 
images of NGC 3256 at its four channels for morpho¬ 
logical and photometric studies. Although Lira et al. 
(2008) (see also Lira et al. 2002) published the IRAC im¬ 
ages and provided nuclear photometries, they did not use 
the 8.0 /xm channel image in their analysis because it is 
slightly saturated at around the N nucleus. The 8.0 ^m 
information is helpful in conjunction with information 
of the 3.6, 4.5, and 5.8 /xm channels to create IRAC 
color-color diagrams for AGN diagnostics (e.g., Lacy et 
al. 2004; Sajina et al. 2005; Stern et al. 2005; Donley et 
al. 2012). Fortunately, the images were taken under the 
HDR (High Dynamic Range) mode in which both short- 
and long-exposure frames are taken within an Astronom¬ 
ical Observation Request (AOR). The short-exposure 
frame is not saturated at the N nucleus. We retrieved 
the Post Basic-Calibrated Data (Post BCD or PBCD) 
of standard mapping observations (AOR: 3896832; PI: 
G. Fazio) from the Spitzer Heritage Archive. This is the 
same data set of Lira et al. (2002, 2008). For the 8.0 /xm 
channel, since short-exposure frames are not recom¬ 
mended for science use in general (IRAC data handling 
book), we confirmed flux calibration consistency between 
the long- and short-exposure frames using photometry of 
field stars. Then the saturated pixels at around the N 
nucleus in the long-exposure frame were substituted with 
corresponding pixels in the short-exposure frame image. 
For the 4.5 /xm channel, additional simple sky subtrac¬ 
tion was made to correct for its tilted background. As 
shown by Lira et al. (2008) for the 3.6, 4.5, and 5.8 /xm 
channels, the two compact nuclei are clearly separated in 
all four channels (Figure 1). 

Aperture photometry was made on each nucleus with 
our IDL program (Table 1). We used an aperture of 278 
diameter for the photometry, and applied aperture cor¬ 
rection following Surace et al. (2004). Flux contribution 
from the host galaxy is estimated within a concentric 
annulus whose inner and outer radii are 270 and 276, 
respectively. The size of the annulus was set as small 
as possible to estimate the flux level near the nucleus 
by minimizing the effect of host galaxy structure. Stan¬ 
dard deviation within the annulus is quadratically added 
to the statistical error of the flux measurement of the 
nuclei. Although our fluxes are systematically slightly 
fainter than those of Lira et al. (2008), the two are con¬ 
sistent because Lira et al. (2008) performed simple aper¬ 
ture photometry with slightly larger apertures than ours 
(376, 376, and 470 diameter apertures for the 3.6 /xm, 
4.5 /xm, and 5.8 /xm images, respectively). We note that 
N nucleus has only about 7%, 7%, 6%, and 7%, and 
the S nucleus has about 4%, 9%, 4%, and 2% in 3.6 /xm, 



An AGN in the NGC 3256 Southern Nucleus 


3 


4.5 /im, 5.8 /im, and 8.0 /im channels, respectively, of the 
flux integrated over the entire galaxy reported by Engel- 
bracht et al. (2008). Flux ratio maps were also made 
for 4.5 /im/3.6 /im, 5.8 /im/4.5 /im, and 8.0 /im/5.8 /tm 
(Figure 2). The point spread function (PSF) size match¬ 
ing was made for these ratio maps using a bright field 
star as a reference. The measured PSF FWHM is about 
2.9 pixel or 1"7. 

2.2. Spitzer IRS 

We processed our own IRS spectra so that the spec¬ 
trophotometry was done in a way consistent with the cor¬ 
responding IRAC photometry for each nucleus, although 
similar spectra had been already published (Brandi et 
al. 2006; Bernard-Salas et al. 2009; Pereira-Santaella et 
al. 2010; Di'az-Santos et al. 2010a; Alonso-Herrero et 
al. 2012). We used the same IRS mapping data set of 
Pereira-Santaella et al. (2010) and Alonso-Herrero et al. 
(2012), and followed almost the same processing meth¬ 
ods as adopted in the earlier studies. The BCD data 
(AOR: 17659904; PI: G. Rieke) of the Short-Low (SL) 
module were retrieved from the Spitzer Heritage Archive. 
Similar mapping data with other modules (Long-Low 
(LL), Long-High (LH), and Short-High (SH)) were not 
used since spatial resolution was worse for the longer 
wavelength data. They are not suitable for our analy¬ 
sis in which separating the two nuclei is essential. The 
standard CUBISM data processing (Smith et al. 2007a; 
Pereira-Santaella et al. 2010) was made, including back¬ 
ground subtraction, hot/rogue pixel identification and 
subtraction, and cube creation. 

Aperture photometry was made with our IDL program 
to extract spectra for each nucleus on the data cubes 
created with CUBISM (Figure 3). We fixed an aper¬ 
ture size of 376 diameter over 5.2-14.5 /rm and applied 
the aperture correction. Flux contribution from the host 
galaxy is estimated within a concentric annulus whose 
inner and outer radii are 270 and 276, respectively. The 
size of the annulus was set as small as possible to esti¬ 
mate flux level near the nucleus by minimizing the ef¬ 
fect of host galaxy structure. However, since the IRS 
spatial resolution (cs 4" FWHM) is larger than that of 
IRAC and is comparable to the distance to another nu¬ 
cleus (—5"), we carefully applied masks to hide another 
nuclear component in the annulus. Standard deviation 
within the annulus is quadratically added to the statis¬ 
tical error of the flux measurement of the nuclei. The 
aperture correction factor was measured with a mapping 
observation of a standard star, HR 7341. The standard 
star data (AOR 19324160) taken in the same observing 
campaign of the NGC 3256 observation were retrieved 
from the Spitzer Heritage Archive, and they were re¬ 
duced by CUBISM with exactly the same parameters. 
We then measured the star with the same photometry 
program and aperture settings, and found the aperture 
correction factor by comparing the measured spectrum 
with the standard-star flux information from the IRS In¬ 
strument Handbook. 

As in the case of the IRAC photometry, those nuclear 
spectrophotometries are about 12-20% of that over the 
entire galaxy (1374 x 13 74; Alonso-Herrero et al. 2012), 
although the overall spectral shapes look similar, i.e., 
all show prominent PAH 6.2, 7.7, 8.6, 11.3 /rm features 
and fine structure lines such as [Ne ii]. We confirmed 


characteristics reported earlier by Pereira-Santaella et al. 
(2010), namely brighter H 2 fluxes, deeper silicate absorp¬ 
tion at 9.7 /xm, and smaller EW of PAH 6.2 /tm at the S 
nucleus (see §3.2 below for our silicate absorption depth 
measurements). 

2.2.1. Flux Matching between the IRAC and IRS 
Spectrophotometries 

The IRS nuclear photometries are about two times 
larger than the IRAC 5.8 /tm and 8.0 /tm nuclear pho¬ 
tometries at both nuclei (Figure 3), and we interpret this 
as a result of more contamination from the host com¬ 
ponent due to larger PSF of IRS. Therefore, in order 
to combine the IRS and IRAC spectrophotometries in a 
consistent way for our subsequent analyses, we adjusted 
the IRAC nuclear photometries for the IRS resolution in 
the following way. For the N nucleus, we simply scaled 
the observed nuclear IRAC fluxes by a factor measured 
at both 5.8 /tm and 8.0 /tm channels by comparing with 
synthetic IRAC fluxes from the IRS spectrum. We uti¬ 
lized a standard photometry tool (spitzer_synthphot.pro) 
from Spitzer Data Analysis Cookbook. Since the IRS SL 
module cannot fully cover the IRAC 5.8 /tm bandpass 
between 4.9 and 5.2 /tm, we extrapolated the spectra 
shortward a little with a help of scaled SWIRE NGC 6090 
SED template (see §3.2 below for more details about the 
template). We found that a single scaling factor on the 
IRAC N nuclear fluxes satisfies matching both 5.8 /tm 
and 8.0 /tm fluxes with the IRS-synthetic photometries, 
and the factor is applied to all four IRAC channel fluxes. 
For the S nucleus, we added scaled IRAC fluxes at the 
N nucleus to the fluxes at the S nucleus to match IRS- 
synthetic photometries at the S nucleus. Here we utilized 
the fact that colors of the N nucleus are very similar to 
those in the surrounding region of the S nucleus (Fig¬ 
ure 2). We calculated the scaling factor by considering, 
again, both 5.8 /xm and 8.0 /xm fluxes. We found that 
a single scaling factor satisfies matching fluxes in both 
channels, and applied the factor to all four IRAC chan¬ 
nels. The adjusted IRAC fluxes for each nucleus are also 
shown in Table 1. Note that the adjusted fluxes are used 
only for analyzing MIR SED together with IRS. 

2.3. T-ReCS Spectrum 

The N-band high-spatial low-spectral resolution spec¬ 
trophotometry of the S nucleus taken by Dfaz-Santos et 
al. (2010a) shows remarkably different signatures just on 
the nucleus, and we therefore include the data in our 
analysis. The spectrum was taken with T-ReCS at the 
Gemini-South Telescope through a 0736-wide slit under 
0735 FWHM seeing condition, providing the best spa¬ 
tial resolution at N band owing to both the diffraction- 
limited 8 m telescope and the seeing-matched narrow slit. 
The spectrum is extracted over 0736 (or four pixel width) 
around the peak at the S nucleus. Dfaz-Santos et al. 
(2010a) performed flux calibration by a standard star 
observation while considering the slit loss effect, and we 
further applied the aperture correction to estimate the 
flux of the central compact source by compensating loss 
of photons outside the extraction aperture box. By as¬ 
suming a Gaussian PSF of the seeing size, we estimated 
and multiplied a factor of 1.30 on the original data pre¬ 
sented by Dfaz-Santos et al. (2010a) in their Figure 3. 
Since error of the spectrum is not explicitly given, we 
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estimated it from a scatter of the data points around the 
fitted model profile (see Figure 3 of Dfaz-Santos et al. 
2010 a). 

2.4. HST NICMOS Images 

In order to extend the shorter wavelength coverage 
of the nuclear spectrophotometries, we also analyzed 
high-resolution HST NICMOS images at NIR. We re¬ 
trieved the combined NIC3 images of FllOW (1.1 /rm 
wide band), F160W (1.6 /rm wide band), and F222M 
(2.2 /xm medium-wide band) from the Hubble Legacy 
Archive. The NIC3 camera provides a larger field-of- 
view and is good for comparison with the IRAC images. 
Their spatial resolutions (PSF size in FWHM) are 0"26 
for all three bands. Flux measurement is based on stan¬ 
dard flux conversion factors for the instrument and the 
filter from the NICMOS Handbook. The HST images 
presented by Lira et al. (2002) and Alonso-Herrero et al. 
( 2002 , 2006a,b) are from a different data set taken with 
the NIC2 camera of the NICMOS instrument. A small 
correction of tilted background was made for each filter 
image by fitting the image with a tilted plane while mask¬ 
ing sources for each filter band. Astrometry is corrected 
with three bright field stars in the 2MASS catalog out¬ 
side of the galaxy, and the global accuracy is about 0 " 2 . 
The F160W and F222M images are shown along with the 
IRAC images in Figure 1. Zoom-up images of all three 
filters around the S nucleus are also shown in Figure 4. A 
compact source marginally appears on the F222M image, 
although it is blended with surrounding structures. This 
source is not clearly identified in FI60W and FIIOW. For 
more NICMOS images for other lines/bands and their de¬ 
scriptions, see Lira et al. (2002) and Alonso-Herrero et 
al. ( 2002 , 2006a). 

3. INFRARED RESULTS 
3.1. Infrared Morphology 

Figure I compares overall morphology of the galaxy in 
all IRAC channels and NICMOS FI60W and F222M im¬ 
ages. In general, IRAC channels (in particular the longer 
5.8 and 8.0 /xm channels) trace star-forming regions and 
the NICMOS bands trace stellar population. The N nu¬ 
cleus is bright and distinct in all IRAC channels. The S 
nucleus is barely visible at 2.2 /xm and is more evident in 
the IRAC channels. The N nucleus is brighter than the S 
nucleus at 3.6, 5.8, and 8.0 /xm, while the two nuclei are 
comparable at 4.5 /xm. Such characteristics were noted 
for 2.2, 3.6, 4.5, and 5.8 /xm by Lira et al. (2002, 2008). 
Kotilainen et al. (1996) also showed that both nuclei are 
equally bright at L’ band (3.6 /xm). 

The flux ratios at the S nucleus are remarkably differ¬ 
ent from the rest due to a compact red source at the S 
nucleus. The S nucleus shows larger 2.2 /xm/1.6 /xm and 
4.5 /xm/3.6 /xm, and smaller 5.8 /xm/4.5 /xm flux ratios 
than those of the N nucleus and most of the host galaxy 
(Figure 2). The large-scale dust lanes known around the 
S nucleus in the optical wavelengths (Zepf et al. 1999; 
Moran et al. 1999; Lfpari et al. 2000; Alonso-Herrero et 
al. 2002, 2006a; Sakamoto et al. 2014) are not visible in 
the MIR flux ratio maps. The red S nucleus stands out 
even more in the HST data when the stellar component 
with complex morphology is subtracted. We modeled the 
stellar component of the host galaxy in the F222M band 


using FI60W and FIIOW images. First we measured 
stellar colors at each pixel by using the two images, and 
we then synthesized the stellar F222M image by linearly 
extrapolating from the FI60W image by considering the 
stellar colors. We scaled the synthesized image by a fac¬ 
tor of 0.74 to match the observed F222M image at well 
outside of the S nucleus. This correction factor is prob¬ 
ably due to color corrections on the filters for very red 
stellar continuum along with the heavy dust lanes around 
the S nucleus (Zepf et al. 1999; Moran et al. 1999; Lfpari 
et al. 2000; Alonso-Herrero et al. 2002, 2006a; Sakamoto 
et al. 2014). We found a distinct compact peak at the S 
nucleus on the F222M image subtracted for its synthetic 
stellar component (Figure 4), indicating an additional 
component at the S nucleus. This compact source coin¬ 
cides with the red peak in the F222M/FI60W image, a 
radio continuum source at the S nucleus (e.g., Norris & 
Forbes 1995; Neff et al. 2003; Sakamoto et al. 2014), and a 
compact (at 0'.'30 FWHM resolution) 8.74 /xm continuum 
source at the S nucleus (Alonso-Herrero et al. 2006b). 

The red compact source at the S nucleus contains an 
unresolved (at 0'.'26 FWHM, or 44 pc) core. We fitted 
the source with one and two circular Gaussian compo¬ 
nents. In the single-Gaussian model, we found a system¬ 
atic radial variation in the residual map in a sense that 
the center has a compact positive residual peak and the 
surrounding region systematically shows negative resid¬ 
uals. In the double-Gaussian model, we assumed that 
the two Gaussian components share the same center po¬ 
sition and that the more compact component is unre¬ 
solved (at 0'.'26 FWHM). The residual map has much less 
systematic structure compared with the single-Gaussian 
model. The reduced are 1.8 and I.O for the single- 
and double-Gaussian models, respectively. Therefore, 
the double-Gaussian model with an unresolved core is 
preferred. The fitting results are summarized in Table 2. 

3.2. MIR Spectra/SEDs and Color-Based AGN 
Diagnostics at Individual Nuclei 

The two nuclei are basically similar to each other in 
the IRS spectra, but we found their noticeable differ¬ 
ences. The IRS spectra of both nuclei show prominent 
PAH and mild silicate absorption features as well as 
some fine structure lines (Figure 3). Such character¬ 
istics are common to LIRGs (e.g., Petrie et al. 2011). 
We compare the observations with the starburst IRS 
template of Brandi et al. (2006), which was generated 
by averaging IRS spectra of local starburst galaxies. 
Since this template is taken with exactly the same in¬ 
strument and module (SL), we can directly compare it 
with our observations. Although the template repro¬ 
duces the observations rather well, we found two no¬ 
table differences between them. One is the depth of 
the 9.7 /xm silicate absorption, which is characterized 
by Asi 9.7 = In (F°b7i^“-‘). Here, is a con¬ 

tinuum flux at 9.7 /xm estimated following Dfaz-Santos 
et al. (2010a). We measured Ssi 9.7 fim = —0.65 ± 0.20, 
— 1.31 ± 0.34, and —0.50 for the N nucleus, S nucleus, 
and the IRS starburst template, respectively. The S nu¬ 
cleus shows significantly deeper absorption than the N 
nucleus, which shows similar depth within uncertainty 
as the starburst template. We note that the N nucleus 
shows slightly smaller Asi 9.7 than the template, al¬ 
though the observed 9.7 /xm flux at the N nucleus is larger 
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than the template scaled at 6.0-8.0 (Figure 3). This 

is because the N nucleus also shows higher continuum at 
>12 fiuY than the scaled template due to systematically 
redder color at > 9 /im, and the estimated continuum 
level at 9.7 /rm is also higher for the N nucleus. Another 
difference is the color at A < 6 ^m. Ratio of our IRS 
observations to the starburst template (Figure 5) shows 
a steep excess at the S nucleus below ~ 6 /im, reaching 
about data-to-template ratio of ~ 2 at ~ 5 /xm. The N 
nucleus does not show such a systematic deviation from 
the template. We also compare the observations with 
a starburst-powered LIRG SED template of NGG 6090 
(log Lir (Lq) = 11.51; Sanders et al. 2003) from Polletta 
et al. (2007) (Figure 3). The template was generated to 
fit the observations with the GRASIL code (Silva et al. 
1998), which is a physical starburst evolution model for 
estimating SEDs. Although this template also roughly 
resembles the spectra of two nuclei, it seems closer to the 
N nucleus at < 6 /im. 

The IRAC SEDs further illustrate the difference be¬ 
tween the two nuclei (Figure 3). At the S nucleus, the 
IRAG flux monotonically increases toward longer wave¬ 
lengths, i.e., no flux drop is found at 4.5 /im. At the N 
nucleus, the IRAG 3.6-4.5 /tm color is almost flat, indi¬ 
cating a contribution from the stellar component, which 
has a blue color at > 1.6 /im (e.g., Sawicki 2002). The 
slope between 4.5 /tm and 8.0 /im is steeper (redder) at 
the N nucleus than at the S nucleus. An excess flux at 

4.5 /im at the S nucleus causes such differences. 

By utilizing IRAC color-color diagrams for AGN di¬ 
agnostics, we found that the S nucleus shows AGN-like 
SED at 3.6-8.0 /im (Eigure 6). Two types of such dia¬ 
grams have been developed for AGN diagnostics. One is 
based on flux ratios of 5.8 /im/3.6 /im vs. 8.0 /im/4.5 /tm 
originally proposed by Lacy et al. (2004). We take color 
boundaries for AGN selection from Donley et al. (2012), 
who modified the flux ratio cut of 5.8 /im/3.6 /tm from 
the original one. Another is based on 5.8 /im-8.0 /tm and 

3.6 /im-4.5 /im colors (in magnitudes) proposed by Stern 
et al. (2005). Although the same IRAC flux information 
is used in both plots, they are not completely consistent 
in identifying photometric AGN candidates. The IRAC 
colors at the S nucleus are in the area for AGNs, which 
include SEDs of pure power-law and hot blackbody of 
300-1000 K. The N nucleus is closer in the IRAC colors to 
starburst-dominated sources (either spiral galaxies, star- 
bursts, or ultra-luminous infrared galaxies (ULIRGs)). 
The difference in colors can be again interpreted as due 
to 4.5 /im excess at the S nucleus. 

3.3. MIR Spectrophotometry Modeling 

It is very likely that the S nucleus is composed of mul¬ 
tiple components with different characteristics. This is 
because the spectrum of the S nucleus at 0'.'36 scale is re¬ 
markably different from the spectra at arcsec resolution 
(Dfaz-Santos et al. 2010a). Within the 0'.'36 aperture, the 
S nucleus shows much deeper silicate absorption without 
detectable PAH. The PAH features are evident in the IRS 
at ^ 3" scale (Dfaz-Santos et al. 2010a; Pereira-Santaella 
et al. 2010, this work), the wider-slit T-ReCS (l'.'3-wide; 
Lira et al. 2008), and the TIMMI2 (l'.'2-wide; Siebenmor- 
gen et al. 2004; Martm-Hernandez et al. 2006) spectra. 
This difference is probably because circumnuclear star¬ 
forming regions dominate the MIR flux within the larger 


apertures. In terms of flux level, all the larger-aperture 
spectrophotometries are comparable to our IRAC 8.0 /im 
nuclear photometry, but the T-ReCS 8.0 /im flux within 
the 0736 aperture is as small as 30% of the IRAC flux. 

We constructed simple composite SED models for the 
S nucleus with and without AGN to examine how addi¬ 
tional AGN contribution can better reproduce the MIR 
observations compared to a pure (but absorbed) star- 
burst model. We fitted the models to both the high 
spatial-resolution T-ReCS spectrum, in particular its 
deep silicate absorption, and the nuclear IRAC and IRS 
spectrophotometries, in particular its prominent PAH 
features. Eor simplicity, only two components are in 
the composites. One of them needs to be compact and 
heavily absorbed to reproduce the very deep silicate ab¬ 
sorption of the T-ReCS spectrum. Another one rep¬ 
resents the extended circumnuclear star-forming region 
which appears to dominate the PAH feature. We made 
two such composite models. The AGN-starburst com¬ 
posite model assumes a heavily absorbed AGN in addi¬ 
tion to surrounding star-forming regions at ~ 3" scale. 
The starburst-starburst composite model assumes heav¬ 
ily and mildly absorbed starburst regions. Emission from 
the heated dusts associated with the compact and heav¬ 
ily absorbed starburst component is not considered in 
the fitting. As a reference, we also constructed a single¬ 
component absorbed starburst model without an AGN. 
Although Alonso-Herrero et al. (2012) already did simi¬ 
lar but more sophisticated SED modeling over the whole 
aperture IRS spectrum, ours is improved in the follow¬ 
ing three points. Firstly, we include the IRAC 4.5 /rm 
flux, below the IRS wavelength coverage, in the fitting 
so that the bluer SED at < 6 /xm of the S nucleus can be 
clearly traced. The IRAC 3.6 /xm flux is not included in 
the fitting because the flux is usually contaminated by 
stellar emission and relative contribution of the stellar 
component with respect to non-stellar component varies 
from one star-forming region to another. Secondly, we 
utilize the nuclear spectrophotometry. Alonso-Herrero 
et al. (2012) used a very large aperture IRS spectrum 
that includes both nuclei and their surroundings. Given 
the only mild differences between the N and S nuclei 
at the IRS spatial resolution (§3.2), it is essential to 
separately analyze the nuclear spectra. Thirdly, we re¬ 
quire our model to fit the high-spatial-resolution T-ReCS 
spectrum while simultaneously fitting the IRAC and IRS 
spectrophotometry. 

We fit the T-ReCS spectrum at 8.0-12.9 /xm (§2.3) 
and the IRAC 4.5 /xm-IRS spectrophotometries at 4.5- 
14.5 /xm (§2.2.1). Eor reference, we also fit only the 
IRS spectra (5.2-14.5 /xm) to evaluate the effect of in¬ 
cluding the IRAC 4.5 /xm photometry. We employed a 
standard minimization technique implemented in the 
MPEIT software package (Markwardt 2009). For the 
starburst component, we adopt the LIRG SED template 
of NGC 6090 from Polletta et al. (2007) since the IRS 
starburst template of Brandi et al. (2006) does not com¬ 
pletely cover the IRAC 4.5 /xm channel. Eor the AGN 
component, we adopt a pure power-law SED (/^ oc xx“) 
with a power-law index (a) of —0.5. The colors/flux ra¬ 
tios of this SED are well within the expected ones for 
AGNs (Eigure 6), and, as we show later, this intrinsi¬ 
cally bluer power-law fits the observed IRAC SED when 
heavily extinct to fit the deep silicate feature of the T- 
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ReCS spectrum. For the extinction curve, we adopt three 
different types and compare the results because, as we 
demonstrate later, a slight difference in silicate absorp¬ 
tion profile can result in quite different best-fit parame¬ 
ters with similarly good statistics for the same T-ReCS 
data. This strongly affects our fit to the IRAC and IRS 
data covering wider wavelengths. We cannot identify 
more appropriate extinction curve without having good 
anchor point (s) of the fit away from the very deep silicate 
absorption. Due to the same problem, we fit the two sets 
of spectrophotometry data one by one because fit to the 
T-ReCS spectrum at 8.0-12.9 /rm hardly constrains our 
fit around 4.5 /rm on the IRAC and IRS spectrophotome¬ 
tries. Therefore, we first fit the T-ReCS spectrum with 
each of the three extinction curves by two SED models 
(either starburst or AGN). During the fit, we fix the red- 
shift, and fit the SED scaling factor and extinction (rep¬ 
resented by optical depth at 9.7 /xm, rg.T^m)- We then 
subtract the fitted model from the IRAC and IRS data, 
and fit the residual spectra with the starburst SED by us¬ 
ing the same extinction curve. We again fix the redshift, 
and fit the SED scaling factor and rg.y^m- We assume a 
screen geometry for both components. We do not con¬ 
sider an effect of radiation from heated dust, although 
such a radiation is expected for a compact dusty region 
surrounding the central energy sources (either starburst 
or AGN). 

We employ the following three extinction curves in the 
fitting. The first one is a theoretical extinction curve of 
Ghiar & Tielens (2006) (hereafter, CT06) for the Galac¬ 
tic center. The second is from a spectral fitting code 
PAHFIT (Smith et al. 2007b). Smith et al. (2007b) in¬ 
troduced a hybrid curve consisting of an observed 9.7 /xm 
silicate absorption profile, an empirical 18 /xm silicate 
absorption profile, and a power-law component. The 
9.7 /xm silicate feature is taken from Kemper et al. (2004) 
based on the observation toward the Galactic center. 
This curve is found to work very well on IRS spectra of 
normal and star-forming galaxies (Smith et al. 2007b; see 
also, e.g., Dfaz-Santos et al. 2010b; Pereira-Santaella et 
al. 2010; Dfaz-Santos et al. 2011). The third one uses the 
same function of PAHFIT but with a different parameter 
set to reproduce the updated extinction curve toward the 
Galactic center by Fritz et al. (2011). The curve of Fritz 
et al. (2011) shows stronger silicate absorption than that 
of CT06, although both curves are for the Galactic cen¬ 
ter. For our fitting purpose, we modified parameters of 
the PAHFIT extinction curve to reproduce the very high- 
resolution extinction curve of Fritz et al. (2011). We set 
the power-law slope of the continuum extinction curve 
zero, and relative strength of the power-law component 
with respect to the 9.7 ^m silicate feature strength (/3; 
Smith et al. 2007b) 0.155. Because the extinction curve 
of Fritz et al. (2011) is almost proportional to that of 
GT06 below the silicate feature, we replaced short side 
of the extinction curve (< 7.3 /xm) with the scaled curve 
of CT06. This gives smaller Av/Tg.y (5.8) than that of 
GT06 (Av/Tg.Y = 9). The three extinction curves are 
compared in Figure 7. Results of the fits are summa¬ 
rized in Figure 8 and Table 3. As we show below, the 
original PAHFIT extinction curve gives much poorer fits, 
therefore the best-fit SED models are not shown with this 
curve in Figure 8. 


3 . 3 . 1 . Absorbed Single-Starburst Model 

In the absorbed single-starburst model (Figure 8 top), 
we modeled the IRAC 4.5 /xm-IRS spectrophotometry 
of the S nucleus only with the absorbed LIRG SED tem¬ 
plate. Although the model reproduces the IRS spectrum 
reasonably well, neither the 4.5 /xm excess nor deep sil¬ 
icate absorption at the 0'.'36 scale can be reproduced. 
Reduced are about 4-6 for all three extinction curves, 
and the IRAC 4.5 /xm flux contributes most to the y^- 

3 . 3 . 2 . AGN-Starburst Composite Model 

In the AGN-starburst composite model (Eigure 8 mid¬ 
dle), we found a good fit, although not statistically satis¬ 
factory, to the T-ReCS spectrum with the AGN template 
by using the extinction curve of CT06 (reduced y^ = 1.4) 
and our modified PAHFIT extinction curve (reduced 
y^ = 1.2). On the other hand, the original PAHFIT 
extinction curve gives poor fit (reduced y^ = 3.6) be¬ 
cause the curve always produces asymmetric extinction 
at both ends of the silicate absorption profile, while the 
observed profile is more symmetric than the curve. We 
note that the best-fit parameters are significantly differ¬ 
ent for different extinction curves (Tg,7=9.4 and 12.7 for 
the extinction curve of GT06 and our modified PAHFIT 
extinction curve, respectively). As explained earlier, this 
is due to slightly different profile of the silicate absorp¬ 
tion feature and absence of fit anchor point (s) outside 
the feature. For the fit to the IRAG 4.5 /xm-IRS spec¬ 
trophotometers, both extinction curve of GT06 (reduced 
y^ = 1.4) and our modified PAHFIT extinction curve 
(reduced y^ = 0.8) resulted in rg.y ~ 0.5-0.7. We note 
that the IRAC 4.5 /xm flux dominates the y^ in a fit 
with the extinction curve of CT06 in a sense that the fit 
under-predicts the observation. 

3 . 3 . 3 . Starburst-Starburst Composite Model 

In the starburst-starburst composite model (Figure 8 
bottom), we found that a heavily absorbed LIRG SED 
also fits the T-ReCS spectrum (reduced y^ = 1.2 for 
the extinction curve of CT06 and our modified PAH¬ 
FIT extinction curve). This is because the silicate ab¬ 
sorption profile of the extinction curve essentially domi¬ 
nates the shape of the absorbed SED within the T-ReCS 
wavelength coverage. Again, the original PAHFIT ex¬ 
tinction curve gives poor fit (reduced y^ = 2.7) because 
of the same reason for the case of the AGN-starburst 
composite model. The fits to the IRAC 4.5 /xm-IRS 
spectrophotometries (reduced y^ ~ 4-6) are much worse 
than those of the AGN-starburst composite model (re¬ 
duced y^ ~ 0.8-1.4), because contribution to the 4.5 /xm 
flux from the heavily absorbed LIRG template is negli¬ 
gibly small. 

3 . 3 . 4 . Summary and Implication of the SED Modeling 

Results 

We found that the AGN-starburst composite model fits 
the MIR observations much better than the starburst- 
starburst composite model since AGN contribution can 
reproduce the 4.5 /xm excess. In particular, the modified 
PAHFIT extinction curve gives better fit to the 4.5 /xm 
flux, and we prefer the results with this extinction curve, 
i.e., Figure 8 right-middle. This model can also roughly 
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explain the NICMOS 2.2 /im flux of the unresolved nu¬ 
clear component, although the flux is not included in the 
fit (Figure 9). In our best-fit model with the modihed 
PAHFIT extinction curve, the AGN contributes ~ 24% 
and ~ 2% of the 6 and 24 /rm fluxes, respectively. As 
for the optical depth on the AGN component, we use 
both results, T9.7 = 9.4 ± 0.4 and 12.7 ± 0.5 from the 
extinction curve of GT06 and our modified PAHFIT ex¬ 
tinction curve, respectively, and adopt the possible range 
of the optical depth of T9.7 = 9-13, since the two extinc¬ 
tion curves provide statistically equally good ht to the 
T-ReCS spectrum. 

We estimate the column density toward the AGN to 
be on the order of A^h ~ 10^^ cm“^. We assumed a 
screen geometry of the dusty region. The optical depth 
at 9.7 jMn over the AGN continuum, T 9.7 = 9-13, corre¬ 
sponds to Ay ~ 80 mag by using the conversion factors 
of Av/t 9.7 = 9 and 5.8 for extinction curves of GT06 and 
our modified PAHFIT, respectively. Further assuming a 
standard conversion factor for the solar neighborhood, 
Nn/ Ay = 1.9 X 10^^ cm“^ mag“^ (Bohlin et al. 1978), 
the neutral Hydrogen column density toward the AGN 
is estimated to be A^h — 1-5 x 10 ^^ cm“^. 

We can also estimate the star formation rate (SFR) 
and IR luminosity of the S nucleus. By subtracting the 
synthetic 8.0 /rm flux of the AGN component (16 mJy), 
we estimate 8 /rm flux from the starburst component of 
our best AGN-starburst composite model is 35 mJy. We 
then estimate SFR to be 0.43 Mq yr“^ by using calibra¬ 
tion of Zhu et al. (2008), which is derived for young (10- 
100 Myr continuous burst) starburst with Salpeter initial 
mass function (0.1-100 Mq) based on IR (8-1000 /rm) 
luminosity-SFR conversion and monochromatic 8 /rm 
luminosity-IR luminosity correlation. This SFR corre¬ 
sponds to Lir, ~ 2.5 X 10® Lq. Zhu et al. (2008) noted 
that correlation between MIR and IR luminosities are es¬ 
sentially the same for AGN-starburst composite galaxies, 
AGNs, and star-forming galaxies. Therefore, the IR lu¬ 
minosities of the AGN and the S nucleus, including the 
AGN and the circumnuclear starburst component, are 
estimated to be ~ 1.1 x 10® Lq and 3.6 x 10® Lq, re¬ 
spectively®. We note that the AGN luminosity is highly 
uncertain because the correlation between the MIR and 
IR luminosities is established for field AGNs found in a 
blank-sky survey (Zhu et al. 2008), and it is not clear 
if the same correlation applies to the heavily absorbed 
AGN at the S nucleus. For comparison, the IR luminos¬ 
ity of the N nucleus estimated in the same way from the 
observed 8.0 ^m flux is cs 1.5 x 10^® Lq. The total in¬ 
frared luminosity at 3-1100 /im (TIR) of the N nucleus is 
also estimated from the whole-aperture TIR luminosity 
by Engelbracht et al. (2008), who used both IRAG and 
MIPS (Rieke et al. 2004) 24 /tm, 70 /tm, and 160 /im 
fluxes. Since the IRAG colors of the N nucleus are very 
similar to that of the rest of the galaxy, and the S nu¬ 
cleus that shows distinct colors consists of only a few 
percent of the whole aperture flux (§ 2 . 1 ), we scale the 
whole-aperture TIR luminosity according to the nuclear 
fluxes at 8.0 /tm. We then found Ltir =(2.0-2.4) xlO^® 

^ Sakamoto et al. (2014) referred to this paper in its submitted 
form for the FIR luminosities. The quoted luminosities are 2.9 x 
10^*^ Lq, ~ 1.5 X 10^*^ Lq, and ~ 5.0 x 10® Lq for the N nucleus, 
the S nucleus, and the AGN in the S nucleus, respectively. During 
the revision of this paper, we revised the numbers as shown here. 


Lq for the N nucleus. 

4. X-RAY DATA 

We analyzed archival data of two Chandra observa¬ 
tions of NGC 3256 to search for X-ray indications of the 
presence of an AGN in the S nucleus. The superb spa¬ 
tial resolution of Chandra enables us to extract X-ray 
spectrum of the S nucleus separated from other emission 
components in the host galaxy. The two Chandra obser¬ 
vations are summarized in Table 4. The observations per¬ 
formed on 2000 Jan 5 and 2003 May 23 are referred to as 
the first and second observations, respectively, hereafter. 
NGC 3256 was located on the back illuminated CCD chip 
ACIS-S3 in both observations. The Chandra Interactive 
Analysis of Observations (CIAO) software package ver¬ 
sion 4.6 combined with the latest calibration database 
(CALDB) version 4.6.3 was used to analyze the Chan¬ 
dra data. The data were reprocessed to generate level=2 
event files using the latest calibration. 

We made a light curve for a source free region in the 
same CCD chip to examine the stability of background 
and discarded time intervals showing high background. 
The background was stable in the second observation, 
while the background rates were relatively high in the 
first. The resulting exposure times after discarding high 
background intervals are 16.2 and 27.2 ksec for the first 
and second observations, respectively, totaling ~ 2.7 
times more exposure time than that in Lira et al. (2002). 
X-ray spectra of the S nucleus were extracted using a 
circular region with a radius of 2.2 pixels or 1"1. Back¬ 
ground spectra were made using a circular region with a 
radius of 5"2 (first observation) or 5'.'9 (second observa¬ 
tion) that does not contain noticeable point sources near 
NGC 3256 and were subtracted from the source spec¬ 
tra. The net counts after the background subtraction are 
shown in Table 4 and the background subtracted spectra 
are presented in Figure 10. The spectra of the S nucleus 
were binned so that each bin contains at least one count. 

Spectral fits were performed by using XSPEC version 
12.8.2. We applied a maximum likelihood method to 
fit the spectra using the modified version of C statistic 
(Cash 1979), in which a Poisson distribution is assumed 
for numbers of counts in each bin. While the absolute 
value of C does not provide goodness of fits, AC can be 
used to examine relative goodness or to generate confi¬ 
dence intervals. The errors are at the 90% confidence 
range for one parameter of interest (AC = 2.7). The 
Galactic absorption Ar = 9.35 x 10®® cm“® (Kalberla et 
al. 2005) calculated by using the tool nh in FTOOLS was 
applied to all the models examined below. The phabs 
model in XSPEC was used for photoelectric absorption, 
in which the cross sections from Balucinska-Church & 
McCammon (1992) with a He cross section in Yan et 
al. (1998) are used. All the spectral components except 
for the Galactic absorption were assumed to be emit¬ 
ted/absorbed at the source redshift. 

5. X-RAY RESULTS 

We fitted the spectra from the two observations simul¬ 
taneously. Since the photon statistics are limited, we as¬ 
sumed same spectral parameters for the two spectra. A 
constant factor was multiplied to the model to represent 
variability between the two observations. This factor was 
fixed at unity for the first observation and left free for 
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the second observation. 

5.1. Power-Law Models 

A simple power law modified by absorption intrinsic to 
the source (model A) was examined first. The resulting 
parameters are summarized in Table 5, and the photon 
index and absorption column density of —0.421 q 27 and 
< 9.7 X 10^^ cm“^ (90% confidence upper limit), respec¬ 
tively, were obtained. Although this model represents the 
overall shape of the spectra, the best-fit photon index is 
extremely flat. If the photon index is fixed at a steeper 
value 1.8, which is typically observed in Seyfert galax¬ 
ies (e.g., Dadina 2008), the C statistic becomes worse by 
AC = 17.7 (Model B). 

An apparently flat spectrum is often explained by a re¬ 
flection dominated spectrum or a combination of heavily 
absorbed and lightly absorbed continua, and these pos¬ 
sibilities are tested below. The former case is expected if 
the primary X-ray emission is obscured by a large amount 
of matter exceeding a column density of ~ 10^^ cm“^ and 
if emission scattered by medium surrounding the X-ray 
source dominates observed spectra. In oder to represent 
a spectrum dominated by reflection by cold matter, the 
pexrav model (Magdziarz & Zdziarski 1995) in XSPEC 
modified by intrinsic absorption was used. The incident 
spectrum was assumed to be power law with an expo¬ 
nential high energy cutoff. The photon index and the 
cutoff energy were fixed at 1.8 and 300 keV, respectively, 
since these parameters were not well constrained from 
the observed spectra. A reflector of a slab shape with 
an inclination angle of 60° was assumed, where 0° cor¬ 
responds to face on. The abundance of the reflector was 
assumed to be solar, where the abundance table of An¬ 
ders & Grevesse (1989) was used. The reflection normal¬ 
ization factor (“ref_refl” parameter) was fixed at —1 so 
that only the reflected continuum emission is obtained. 
This model (Model C) represents the observed spectral 
shape, and the best-fit parameters are shown in Table 5, 
though the absence of Fe-K emission line is not compat¬ 
ible with this refection dominated model as discussed in 
§ 6 . 1 . 2 . 

Another possibility to explain the very flat continuum 
is a multicomponent model consisting of heavily and 
lightly absorbed continua. We assumed a power law con¬ 
tinuum with a fixed photon index of 1.8 as an incident 
spectrum and the absorption column densities of Anp 
and Ah ,2 for light and heavy absorbers, respectively. 
This model is expressed as [yg-'T(£;)JVH ,2 _|_ 

(1 — /)] A E~^^ where A, /, A, and F are photon 

energy, photoelectric absorption cross section, fraction of 
continuum emission absorbed by Ah. 2 , normalization of 
power law, and photon index, respectively (Model D). 
The observed spectra are well represented by this model. 
The result of this fit is also in Table 5. The best-fit model 
is shown in Figures 10 and 11. The measured column 
density of the heavy absorber (7^3® x 10^^ cm“^) is con¬ 
sistent with the MIR measurement of an AGN absorption 
from the AGN-starburst composite model (~ 1.5 x 10^^ 
cm“^) (§3.3.4). 

The constant factors multiplied to the models obtained 
from the fits are summarized in Table 5. The best-fit val¬ 
ues for the models examined above are in the range from 
0.92 to 0.95, and the error ranges contain unity. There¬ 


fore, flux variability between the two observations is not 
significant. Observed fluxes and luminosities corrected 
for absorption in the 2-10 keV band are also shown in 
Table 5 except for the model B, for which C statistic was 
significantly worse than those for other models. 

The observed spectra do not show an indication of an 
Fe-K emission line. We calculated an upper limit on the 
EW of a fluorescent Fe-K line at 6.4 keV by adding a 
Gaussian spectral component to the models examined 
above. The center energy was fixed at 6.4 keV. The line 
width CT is assumed to be 10 eV, since previous observa¬ 
tion of type 2 AGNs show that Fe-K line width is much 
narrower than the energy resolution of ACIS (Shu et al. 
2011). The upper limit on the line EW is shown in Ta¬ 
ble 5. 

5.2. Thermal Emission Models 

We also examined a thermal emission model. The 
APEC model (Smith et al. 2001) was used as emission 
from optically-thin plasma in collisional ionization equi¬ 
librium. A combination of multiple components, one of 
which is heavily absorbed, is required to explain the ap¬ 
parent flat spectrum. We therefore examined a model 
consisting of unabsorbed and absorbed APEC compo¬ 
nents represented by the expression apec(Ti) 
apec(T 2 ), where Ti and T 2 are temperatures for unab¬ 
sorbed and absorbed APEC components, respectively 
(Model E). Only lower bounds for allowed temperature 
ranges for these two components were obtained as 1.2 
keV and 9.4 keV, respectively. The intrinsic column 
density for the absorbed component was obtained to be 
Ah = 5.1f'^7 g X 10^^ cm“^. This model resulted in the 
value of C statistic 77.6 for 88 degrees of freedom and 
describes the shape of the observed spectra, where the 
plasma temperatures kTi and kT 2 were pegged at 64 keV. 
X-ray spectra of starburst galaxies generally show emis¬ 
sion from thermal plasma with a temperature kT = 1 
keV or less (Ptak et al. 1999; Strickland et al. 2004; 
Tiillmann et al. 2006), and the high temperatures ob¬ 
tained from the fits are unusual. We therefore exam¬ 
ined whether the presence of low temperature (kT < 1 
keV) plasma is consistent with the observed spectrum 
by multiplying extra absorption to the model examined 
above. This trial model is expressed by [ 

apec(Ti) -p apec(r 2 )], and we found that kTi 

of 1 keV or less is allowed for the component representing 
the low energy part of the spectra (Model F). If kTi= 
1 keV is assumed, the best-fit Ah,i is 1.3 x 10^^ cm“^, 
and C=77.3 is obtained for 88 degrees of freedom. The 
resulting spectral parameters, observed fluxes, and ab¬ 
sorption corrected luminosities for the thermal models 
are summarized in Table 6. 

6. DISCUSSION 

6.1. Evidence of AGN in the S Nucleus 
6.1.1. In Infrared 

The IRAC colors at the S nucleus indicate a non- 
starburst-like SED (§3.2). Specifically, the excess of the 
4.5 /xm flux makes 8.0 /xm/4.5 /xm flux ratio smaller and 
[3.6] —[4.5] color larger in both diagnostic diagrams in 
Figure 6, bringing the S nucleus away from starburst- 
dominated galaxies in the observed colors. Such col- 
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ors are found only at the S nucleus in the IRAC im¬ 
ages (§3.1). The MIR spectrum of the flux ratio between 
the two nuclei (Figure 5a) can be used to examine pres¬ 
ence of an AGN in a way least dependent on the SED 
templates. It shows a complicated trend of both the ex- 
cessed 4-6 fj,m flux and the deeper 9.7 /im absorption 
at the S nucleus. The ratio changes in a way different 
from the dust extinction curves (Figure 7) and, therefore, 
the two nuclei must have different intrinsic SED shapes. 
This conclusion is independent of the detailed shape of 
the extinction curve. Such an excess at 4-6 fini cannot 
be reproduced only from starburst-dominated SEDs be¬ 
cause the observed 8.0 /xm/4.5 /im flux ratio is smaller 
than, and the observed [3.6] —[4.5] color is larger (red¬ 
der) than, those of starburst-dominated galaxies (Ruiz 
et al. 2010). We demonstrated with our AGN-starburst 
composite SED model that the AGN contribution can 
reproduce the observed 4.5 /im enhancement (§3.3.4). 

We discovered a very compact 2.2 /im source at the S 
nucleus, and its unresolved core (at 0'.'26 FWHM) seems 
a NIR counterpart of the compact (at O'.'SO FWHM res¬ 
olution) MIR core (Alonso-Herrero et al. 2006b). The 
2.2 /im core coincides with the compact MIR and radio 
sources at the S nucleus (§3.1). Both the 2.2 /im core 
and the MIR core show a comparable source size. Al¬ 
though this 2.2 /im component is not included in our 
SED fitting together with the T-ReCS spectrophotome¬ 
tries, the heavily extinct power-law AGN SED that is 
fitted to the MIR data roughly predicts the 2.2 /im flux 
(§3.3.4). Therefore, it is very likely that the same power- 
law component, or hot blackbody-dominated component 
(see below), dominates the NIR-MIR SED of the nucleus 
at < 50 pc scale. 

Although our SED model fitting prefers a model with 
AGN, the nuclear IRAC colors can be also represented 
by a blackbody of ~ 600 K (Eigure 6). Although we 
neglected dust emission in our SED modeling of the 
starburst-starburst composite model (§3.3), such dust 
emission might dominate the SED around 4.5 /im to ac¬ 
count for the 4.5 /xm excess. Verley et al. (2007) stud¬ 
ied SEDs of embedded star-forming regions by consid¬ 
ering hot dusty shell with radiative-transfer calculation 
with DUSTY code (Ivezic & Elitzur 1997). No PAH 
emission is considered within DUSTY. The young star 
clusters are modeled with the starburst99 (Leitherer et 
al. 1999) stellar-population model. Among the parame¬ 
ter space they explored, which was designed to cover a 
range of typical extragalactic H ii regions, the S nucleus 
is closer in IRAC colors to cases of temperature of the 
inner surface of the dusty shell Tin = 600 K and total 
visual extinction Ay = 46.4 mag. A heavy 9.7 /xm sil¬ 
icate absorption is also reproduced when the extinction 
is very high (Ay ~ 100 mag). In such configurations of 
the dusty shell, large extinction suppresses stellar SED 
at 3.6 /im, and the color becomes similar to a pure black- 
body of T = 600 K. Luminosity of the 600 K blackbody 
component is ~ 3 x 10® Lq if the blackbody dominates 
the observed nuclear IRAC SED. If we assume a shell 
of dusts emitting such blackbody, radius of the shell is 
< O.I pc. 

MIR SEDs with dominating hot blackbody are un¬ 
usual for star-bursting galactic nuclei. Some types of 
star-forming objects are known to show SEDs similar to 
the IRAC SED at the S nucleus, but they are likely not 


responsible for the S nucleus for the following reasons. 
Hot blackbody-dominated SEDs are found in blue com¬ 
pact dwarf galaxies (Hunt et al. 2005) and ultra-compact 
H II regions (e.g., Churchwell 1990, 2002). In blue com¬ 
pact dwarf galaxies, vigorous star formation is clearly 
seen in optical wavelength, and the MIR SEDs resemble 
a blackbody because of lack of prominent PAH emission 
due to low metallicity. Because the S nucleus is most 
likely a nucleus of a massive disk galaxy under merging 
process (e.g., English et al. 2003; Sakamoto et al. 2014) 
and the metallicity of the nuclear region is about solar 
(e.g., Storchi-Bergmann et al. 1995; Lipari et al. 2000), 
the same explanation for blue compact dwarf galaxies un¬ 
likely applies to the S nucleus. In addition, typical X-ray 
luminosity of such galaxies is about a few 10®® erg s“® 
(Kaaret et al. 2011; Thuan et al. 2014), which is about an 
order of magnitude fainter than the observed luminosity 
of the S nucleus. Ultra-compact H ii regions are manifes¬ 
tations of newly formed massive stars that are still em¬ 
bedded in their natal molecular cloud, and hence dusty 
cocoon (Ghurchwell 1990). It is compact (<I pc) and 
emits mostly at infrared from hot dust heated by central 
O star. Strong silicate 9.7 /xm absorption is sometimes 
seen. Although they are among the most luminous EIR 
sources in the Galaxy and their MIR spectral character¬ 
istics are similar to those of the S nucleus, even the most 
luminous ultra-compact H ii regions in a classical sample 
of Wood & Churchwell (1989) are ^ 4 x 10® times less 
luminous than that of the S nucleus at ~ 10 /xm. Their 
X-ray luminosities (typically a few 10®® erg s“®; Tsuji- 
moto et al. 2006) are also much fainter than the observed 
luminosity of the S nucleus. Nuclear starburst galax¬ 
ies and (U)LIRGs hardly show such hot dust component 
(Marshall et al. 2007; Armus et al. 2007; da Cunha et al. 
2010). In a sample of representative local star-forming 
ULIRGs, da Cunha et al. (2010) requires no component 
hotter than 250 K in their UV-FIR SED fitting. Armus et 
al. (2007) pointed out from their detailed analysis of the 
infrared (I-IOOO /xm) SEDs for local ULIRGs with and 
without AGNs that detection of hot dust at T > 300 K 
in a nuclear spectrum provides indirect evidence for a 
buried AGN. 

In summary, our infrared analysis strongly prefers pres¬ 
ence of an AGN at the S nucleus. In order to better con¬ 
strain the nuclear activities, sophisticated SED decompo¬ 
sition as Marshall et al. (2007) did or physical SED mod¬ 
eling with dust radiation-transfer consideration would be 
necessary. Additional high spatial-resolution photometry 
data to fill a wavelength gap between NICMOS and N 
band as well as spatially-resolved photometry of the S 
nucleus at longer infrared wavelength would greatly help 
such modeling and analysis. 

6.1.2. In X-Ray 

Chandra imaging of NGC 3256 shows an X-ray source 
at the position of the S nucleus. Its spatial extent is con¬ 
sistent with the PSE of Chandra (Lira et al. 2002). We 
first discuss possible origins of the X-ray emission assum¬ 
ing that X-rays are coming from one point source. The 
X-ray spectra are extremely hard; if a simple power law 
model modified by intrinsic absorption is applied, the 
best-fit photon index of —0.43 is obtained (Model A). 
Such a flat spectrum is unusual for primary X-ray emis¬ 
sion and suggests a contribution of reprocessed and/or 
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absorbed emission. If the primary source is hidden be¬ 
hind optically thick matter with a column density greater 
than Nn ~ 10^'^ cm“^ and the observed spectrum is dom¬ 
inated by emission scattered by cold matter, the resulting 
spectrum becomes very hard. Such a spectrum is referred 
to as reflection dominated. This situation is observed in 
Compton-thick AGNs, for which an absorption column 
density exceeds 1.5 x lO^"* cm“^ (e-g-, Comastri 2004). 
In this case, a strong fluorescent Fe-K emission line with 
an EW greater than ~ 700 eV is seen (e.g., Guainazzi 
et al. 2005; Fukazawa et al. 2011). The observed spec¬ 
tra, however, do not show an indication of Fe-K emission 
line, and the upper limit on the EW is 190 eV (for model 
C). This limit is inconsistent with a reflection dominated 
spectrum and the presence of a Compton-thick nucleus 
is ruled out. 

A dual absorber model (Model D) also provided a good 
description of the spectrum. Such spectra are seen in 
AGNs obscured by Compton-thin matter. The fraction, 
/ = 0.94, of the continuum absorbed by a large col¬ 
umn density (7 x 10^^ cm“^) is in the range typically 
observed in Seyfert 2s (Turner et al. 1997; Noguchi et al. 
2010). The less absorbed emission component is often 
interpreted as emission scattered by ionized medium in 
the opening part of the putative obscuring torus (e.g.. 
Turner et al. 1997). The upper limit on the EW of 
an Ee-K emission line (550 eV) is consistent with that 
expected for the absorption column density of 7 x 10^^ 
cm“^ (EW~ 50-150 eV; Turner et al. 1997; Guainazzi 
et al. 2005; Eukazawa et al. 2011). The luminosity in the 
2-10 keV band corrected for absorption is estimated to 
be 1.5 X 10'^° erg s“^. The error on the absorption col¬ 
umn density introduces uncertainties in the correction of 
absorption, and the allowed range for the 2-10 keV lu¬ 
minosity is (1.2-2.9)xl0‘*° erg s“^. This luminosity is in 
the range for low-luminosity AGNs (e.g.. Ho et al. 2001; 
Terashima et al. 2002; Cappi et al. 2006). Thus the ob¬ 
served properties are in accordance with a low-luminosity 
AGN obscured by Compton-thin matter. 

The luminosity of the X-ray source at the S nucleus is 
much higher than that of a single X-ray binary containing 
a stellar mass black hole or a neutron star. If the X-ray 
source is a single compact object, its mass should be 
greater than 
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(Gonzalez-Martm & Vaughan 2012), and the absence of 
variability is consistent with X-ray emission being from 
a low-luminosity AGN. On the other hand, the absence 
of variability might be explained by the possibility that 
the X-ray emission comes from multiple sources. If there 
are tens of stellar mass black holes well within the size 
of the PSF of Chandra (0749 or 83 pc at the distance 
of NGC 3256), the luminosity measured could be ex¬ 
plained. The very hard spectra observed, however, are 
not compatible with any spectral states of stellar mass 
black hole binaries (McClintock & Remillard 2006; Done 
et al. 2007) and superposition of X-ray spectra suffering 
from different absorption column densities is required. A 
large number of X-ray binaries are expected in starburst 
galaxies. The integrated luminosity of high mass X-ray 
binaries, which have a hard X-ray spectrum, formed by 
starburst activity is related to SER (Ranalli et al. 2003; 
Grimm et al. 2003; Gilfanov et al. 2004). The scaling 
law of Grimm et al. (2003) between SER and integrated 
X-ray luminosity (L 2-10 kev) for low SER values (< 4.5 
M 0 yr“4) predicts T 2-10 keV of 6 x 10^® erg s”^ for SER 
of ~ 0.43 M 0 yr“4 (§3.3.4). The observed luminosity 
is about a factor of ~ 20 larger than this expectation. 
Chandra images of nearby starburst galaxies indeed show 
many point sources (e.g., Griffiths et al. 2000; Strickland 
et al. 2000; Bauer et al. 2001). Such point sources, how¬ 
ever, are distributed over a region of several hundreds 
of parsecs, which is much larger than the limit on the 
source size of the S nucleus. A single nuclear source is 
therefore more conceivable, although the possibility of 
superposition of multiple sources cannot be completely 
excluded. 

A two component thermal plasma model also explains 
the shape of the observed spectra (§5.2). Thermal plasma 
emission with a temperature lower than ~ 1 keV is com¬ 
monly observed in starburst galaxies (e.g., Strickland et 
al. 2004). Such a component is extended to the scale of 
host galaxy and intrinsic absorption is generally small. 
Our spectral fit indicates that emission from plasma with 
kT < 1 keV should be significantly absorbed by a column 
density of 1.3 x 10^^ cm“^ or higher if such a component 
exists. Thus the observed properties are unusual to be 
interpreted as starburst galaxies in terms of the small 
source size and \ow-kT plasma conhned and hidden be¬ 
hind absorbing matter. 

6.1.3. Arguments at Other Wavelengths 


where Asdd, '«2-io keV, L 2-10 keV are the Eddington ra¬ 
tio, bolometric correction factor, and intrinsic luminosity 
in 2-10 keV, respectively. Note that we assumed an Ed¬ 
dington luminosity of 1.5 x 10^® erg s“4 for one solar 
mass, which is based on an assumption of a Hydrogen- 
to-Helium ratio of 0.76:0.23 by weight. The bolometric 
correction factor «: 2 -iokeV is known to be dependent on 
a luminosity and an Eddington ratio. Studies of AGN 
SEDs show that K 2 -iokeV is 10-20 for low-luminosity 
AGNs accreting at an Eddington ratio of <0.1 (Vasude- 
van & Fabian 2009; Ho 2008). 

The X-ray fluxes of the S nucleus for the two observa¬ 
tions with a 3 year interval are consistent with each other 
within errors. A systematic analysis of AGN X-ray light 
curves of a large sample shows that a high percentage of 
low-luminosity AGN do not show significant variability 


In addition to our MIR and X-ray evidence for the 
AGN, recent radio and NIR observations also suggest 
the AGN activity in the S nucleus. Sakamoto et al. 
(2014) found a highly collimated bipolar jet-like outflow 
of molecular gas from the S nucleus using ALMA. It ex¬ 
tends up to 4" (700 pc) from the nucleus and is asso¬ 
ciated with a bipolar spur of radio continuum emission. 
They inferred from the morphology and kinematics of 
the molecular outflow and from the radio feature a jet¬ 
driving AGN in the S nucleus. They also found that the 
continuum spectral slope at 860 ^m is flatter at the S nu¬ 
cleus than at the N nucleus. More synchrotron emission 
from an AGN in the S nucleus would explain the differ¬ 
ence. Similar bipolar outflow from the S nucleus was also 
detected in the NIR H 2 1-0 S')!) line by Emonts et al. 
(2014), and was suggested to be AGN-driven from the 
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energetics and high-mass loading factor. 

In addition to the X-ray to MIR SED (Lira et al. 
2002; see also §6.3 below), the X-ray to radio SED has 
been used to distinguish types of nuclear activities be¬ 
cause these wavelengths are least sensitive to extinction 
and contamination by stellar activities is minimal there 
(Terashima & Wilson 2003). Neff et al. (2003) found the 
ratio of 6 cm radio to 2-10 keV X-ray luminosities, i?x 
(= Lr/Lx = vL, (5 GHz)/L 2 -io kev)= 1 X 10-2, for 
the compact source in the S nucleus from their VLA ob¬ 
servations and the Chandra X-ray results of Lira et al. 
(2002). By comparing with the ratios of various kinds 
of Galactic and extragalactic sources, they argued that 
the compact source is most likely a low-luminosity AGN, 
although a possibility of a collection of supernova rem¬ 
nants cannot be rejected. Our new X-ray luminosity is 
only 26% larger than that of Lira et al. (2002), and the 
revised i?x remains in the range of low-luminosity AGNs 
as Neff et al. (2003) discussed. 

6 . 2 . Distribution of Dusty Material in and around the S 

Nucleus 

Our absorption magnitude and column density toward 
the AGN in the S nucleus, Ay ~ 80 mag and A^h ~ 
cm-2, are much larger than previous estimates at NIR- 
MIR wavelengths. They are Ay = 5.3 and 10.7 mag 
from the JHK'L' continuum and line photometries, re¬ 
spectively (Kotilainen et al. 1996), 16 mag from the HST 
NICMOS H~K nuclear color measurement (Lira et al. 
2002), and 15 ± 5 mag from the HST NIGMOS line flux 
ratio (Alonso-Herrero et al. 2006a). Also, Lira et al. 
(2008) estimated Ay ~ 10 mag on the basis of MIR SED 
fitting within 3'.'6-4'.'0 apertures. It seems that the pre¬ 
vious measurements of Ay = 5-15 mag (i.e, A^h ~ 
cm- 2 ) correspond to the extended dust lanes over 1 kpc 
scale (e.g., Zepf et al. 1999; Moran et al. 1999; Lfpari 
et al. 2000; Alonso-Herrero et al. 2002, 2006a; Sakamoto 
et al. 2014) or dusty circumnuclear star-forming region 
(§3.3). In our MIR SED analysis, we applied mild ex¬ 
tinction corresponding to Ay ~ 5-8 mag on the star- 
burst component of the AGN-starburst composite model 
(§3.3.2). This component is most likely spatially ex¬ 
tended because Alonso-Herrero et al. (2006b) reported 
a slightly extended structure in a 0'.'30 FWHM resolu¬ 
tion (51 pc) image of the S nucleus at 8.74 /xm. We 
also found a slightly extended component, besides an 
unresolved core, in the NIGMOS 2.2 /xm image (0'.'75 
FWHM or 128 pc) (§3.1). In our X-ray analysis, we 
showed that the lightly absorbed component of the dual¬ 
component model suffers from A^h,i = x lO^^ 

cm-2 (§5.1), matching the values from previous and our 
measurements at NIR-MIR. Additional heavy extinction 
on the AGN is required to account for both the deep sil¬ 
icate absorption at the 0736 scale and the heavy absorp¬ 
tion on another component of the dual-component model 
in the X-ray analysis. Such a heavily absorbed region is 
most likely within the central < 0.5" of the S nucleus. 
At MIR, the silicate absorption becomes much shallower 
at arcsec scale (Dfaz-Santos et al. 2010a; §3.3). At X- 
ray, high spatial-resolution Chandra spectra show that 
the column density on a heavily absorbed component 
(A^h ,2 = X cm- 2 ) is much larger than that on 
the power-law component in the whole-aperture XMM- 


Newton spectrum (0.14 ± 0.03 x 10^^ cm“2; Pereira- 
Santaella et al. 2011). For comparison, from high spatial- 
resolution ALMA molecular-gas imaging, Sakamoto et 
al. (2014) obtained a column density of A^h = lO^^-lO^^ 
cm“2 to the S nucleus in their 075 (80 pc) beam. 

Both MIR and X-ray analyses strongly suggest that 
large amount of dusts distributes immediately around 
the AGN in a form of thick and smooth shell or 
torus. ULIRGs generally show deeper silicate absorption 
('S'g.T^m 'C —1) than in AGNs (e.g., Hao et al. 2007), and 
such deep absorption requires that the energy source is 
deeply embedded in dust that is both optically and ge¬ 
ometrically (along radial direction) thick, as well as ge¬ 
ometrically smooth (Levenson et al. 2007; Nenkova et 
al. 2008a,b; Sirocky et al. 2008). For example, slab ge¬ 
ometry cannot make the absorption S'g.Ty^m < — l.I even 
with Ty = 1000, because one needs temperature gradi¬ 
ent within the dusty region and the illuminated surface 
of the dust should be hidden from our direct views (Lev¬ 
enson et al. 2007). Also, clumpy medium cannot make 
the absorption very deep because illuminated surface of 
individual clumps can be directly seen/illuminate dark 
surface of other clumps (Nenkova et al. 2008a,b). The 
observed upper limit of the silicate absorption depth at 
the 0736 resolution, S'g.y^m < —3 (Dfaz-Santos et al. 
2010 a), where the bottom of the feature was not deter¬ 
mined due to very strong absorption and the sensitiv¬ 
ity, is already in the range of the most heavily absorbed 
sources in the local universe (Hao et al. 2007). Our MIR 
SED analysis implies a much deeper depth (S'g.y^m = —9 

-13 where S'g.T/^m = —Tg .7 in our assumed geometry). 

This strengthens a need for optically and geometrically 
thick dusty region immediately surrounding the AGN. 
On the other hand, our X-ray spectrum analysis with a 
dual-component model suggests that a classical picture of 
type-2 AGNs with a dusty torus immediately surround¬ 
ing a supermassive black hole (e.g., Antonucci 1993) ap¬ 
plies to the S nucleus. Firstly, the elevated column den¬ 
sity at the S nucleus is compatible with a typical dusty 
torus of Compton-thin type-2 AGNs (e.g., Awaki et al. 
1991; Turner et al. 1997; Guainazzi et al. 2005; Fukazawa 
et al. 2011). Secondly, the model also suggests that ~ 6 % 
of the primary X-ray spectrum is seen unabsorbed by the 
heavy absorber, and this component can be interpreted 
as emission scattered by ionized medium in the opening 
part of the putative obscuring torus around the AGN. 
Such phenomenon is commonly seen in typical Seyfert 2s 
(Turner et al. 1997; Noguchi et al. 2010). To reproduce 
deep silicate absorption with such a torus, the torus must 
be geometrically thick enough along the height direction, 
as well as along the radial direction, or our viewing angle 
must be close enough to the equatorial plane of the torus, 
so that the inner illuminated surface is hidden from our 
line-of-sight. 

6.3. Revisiting X-Ray-to-MIR SED 

The S nucleus is known to be X-ray-quiet, and the 
fact has been considered to be evidence against an AGN 
there (Lira et al. 2002), although the galaxy as a whole 
is among the most luminous X-ray sources without con¬ 
firmed AGN in the local universe (Moran et al. 1999; 
Lira et al. 2002; Pereira-Santaella et al. 2011). Lira et 
al. (2002) revealed a compact X-ray source at the S nu¬ 
cleus with the first Chandra observation and reported the 
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neutral Hydrogen column density obscuring the nucleus 
to be -/Vh = 5 X 10^^ cm“^ (best fit)-l x 10^^ cm“^ (ac¬ 
ceptable) . By using a radio-FIR-MIR-X-ray SED of the 
S nucleus, Lira et al. (2002) showed that the observed X- 
ray luminosity is higher than that of typical starbursts 
but is at least two orders of magnitude lower than ex¬ 
pected for a classical Seyfert nucleus. Alexander et al. 
(2005) plotted the data of Lira et al. (2002) for the to¬ 
tal aperture in the correlation diagram between IR and 
(absorption-corrected) X-ray luminosities together with 
well-studied local starbursts and AGNs. They demon¬ 
strated that the X-ray luminosity of NGC 3256 between 
0.5 keV and 8 keV is well below the range of AGNs but 
the galaxy is slightly more X-ray luminous than the star- 
burst galaxies. 

Our AGN-starburst composite model indicates that 
the observed SED is heavily contaminated at MIR by 
the circumnuclear star formation, and we found that the 
AGN component alone is consistent with typical AGNs 
in terms of X-Ray-to-MIR SED. Horst et al. (2006, 2008) 
demonstrated that contamination by the circumnuclear 
starburst at MIR must be removed when evaluating AGN 
SEDs. By using high spatial-resolution MIR photome¬ 
tries, Horst et al. (2006, 2008) and Gandhi et al. (2009) 
showed that the ratio of the MIR luminosity at 12.3 /im 
to the absorption-corrected X-ray luminosity at 2-10 keV 
scatters within a range of log L 12.3 /im/L2-io keV = —0.5 
to 1.5, with a mean of 0.6I±0.3I among both type-1 and 
-2 Seyfert nuclei. In the case of the S nucleus, we adopt 
the T-ReCS flux at 12.3 /xm within the 0"36 aperture 
(Diaz-Santos et al. 2010a) and the absorption-corrected 
X-ray luminosity of the dual-component model (Model 
D; §6.1.2). Then we find log L 12.3 /im/L 2 -io keV ~ 0.8- 
1.2 (= 1.1 with the best-fit 2-10 keV luminosity) for the 
AGN in the S nucleus. It is in the range of ratios that 
Horst et al. (2008) and Gandhi et al. (2009) obtained for 
Seyfert nuclei. Therefore, we conclude that the appar¬ 
ent X-ray quietness of the S nucleus should no longer be 
considered as evidence against AGN in the nucleus. 

6.4. Comparison with Previous Results 

Although our conclusion of the presence of an AGN in 
the S nucleus appears contradictory to most of the earlier 
studies, they are consistent with each other if we consider 
different aperture sizes and contamination from circum¬ 
nuclear star formation. Lira et al. (2002, 2008) noted the 
4.5 fj,m excess at the S nucleus among the three IRAG 
channels. Lira et al. (2002) also noted, by adding in¬ 
formation of ground-based JHK'L'N-hand photometry, 
that the S nucleus is in the 1-10 /xm SED in vF^, flatter 
than typical starburst galaxies and is closer to AGNs. 
However, they concluded that the AGN contribution is 
insignificant, if any. This is because in Lira et al. (2002) 
their FIR-MIR-X-ray SED for the S nucleus is not con¬ 
sistent with the SEDs of typical Seyfert 2 nuclei, and in 
Lira et al. (2008) their T-ReCS iV-band spectrum with a 
1"3 slit shows evidence for star formation (e.g., PAH fea¬ 
tures). In our work, we showed that the circumnuclear 
star formation dominates the MIR flux at a ^ 3" scale, 
and it is not surprising that the T-ReCS spectrum from 
the 1"3 slit shows PAH features because the T-ReGS flux 
and the IRS flux within the 3"6 aperture are almost the 
same (§3.3). We also showed that the X-ray-to-MIR lu¬ 
minosity ratio is consistent with typical AGNs if we con¬ 


sider only the AGN contribution at MIR after excluding 
contribution from the circumnuclear star-forming regions 
(§6.3). 

Our deeper Chandra spectrum enabled us to better 
constrain the X-ray characteristics owing to improved 
statistics. On the basis of their first Chandra observa¬ 
tion, Lira et al. (2002) presented spectral analysis of a 
composite X-ray spectrum of three hard sources includ¬ 
ing the S nucleus, and constrained the spectral shape of 
the S nucleus using X-ray hardness. We used two Chan¬ 
dra data sets totaling ~ 2.7 times more effective exposure 
time, and were able to obtain spectra of the S nucleus 
alone. The relatively flat spectra we obtained are quali¬ 
tatively consistent with the results by Lira et al. (2002). 
They assumed an absorbed power law model to constrain 
the allowed ranges of the photon index and absorption 
column density. Our spectral fits indicate that a single¬ 
component absorbed power law model gives an extremely 
flat photon index (F < 0.25), otherwise the quality of fit 
becomes significantly worse. The absence of a strong Fe- 
K fluorescent line in our spectrum strongly prefers the 
model consisting of two components, one of which is ab¬ 
sorbed by Gompton-thin matter (iVn ~ 7 x 10^^ cm“^), 
rather than a single absorbed power law. 

The integrated spectrum of NGC 3256 measured with 
XMM-Newton provided a strong constraint on the flux of 
an Fe-K fluorescent line (Pereira-Santaella et al. 2011). 
Their upper limit on the flux is 3.5 x 10“^® erg cm“^ 
s“^ or 3.4 X 10“^ photons cm“^ s“^. Our limit on the 
Fe-K flux (3.6 x 10“^ photons cm“^ s“^ for Model D) is 
almost the same as that obtained by Pereira-Santaella et 
al. (2011). We obtained spectra of the S nucleus alone, 
in contrast to the integrated spectra previously reported, 
and succeeded to set a limit on an Fe-K EW (< 190 eV 
for Model C and < 550 eV for Model D) that excludes 
the possibility of the presence of a Compton-thick AGN. 

The AGN in the S nucleus is energetically much less im¬ 
portant at MIR than the rest within NGC 3256 (§3.3.4), 
and previous estimates of the fractional AGN contribu¬ 
tion within the galaxy seem not sensitive enough to de¬ 
tect the AGN. Alonso-Herrero et al. (2012) made SED 
modeling of the IRS spectrum extracted over a kpc-scale 
region to explore possible AGN contribution. They con¬ 
cluded that the AGN contributions to the whole-aperture 
MIR (at 6 and 24 /xm) and IR luminosities are very small 
(< 5% and < 1%, respectively). Our estimates of the 
AGN contribution to the S nucleus at ~ 3" scale are 
~ 24% and ~ 2% at 6 and 24 /xm, respectively (§3.3.4). 
Since the entire-galaxy IRS spectrum of Alonso-Herrero 
et al. (2012) is about 8 times brighter than our nu¬ 
clear ones (§2.1, 2.2), the AGN contributes to ^ 3% and 
^ 0.3% of the whole-aperture MIR luminosities at 6 and 
24 /xm, respectively. Therefore, our results are consistent 
with those of Alonso-Herrero et al. (2012). 

The absorption corrected X-ray luminosity of the S nu¬ 
cleus for the dual absorber model with Seyfert-like power- 
law continua (Model D), which we conclude most plau¬ 
sible, is 1.5 X 10"^° erg s“^. It is about 20% of the total 
X-ray luminosity integrated over NGC 3256, 7.4 x 10^° 
erg s“^, measured with XMM-Newton (Pereira-Santaella 
et al. 2011). 

7. SUMMARY AND CONCLUSIONS 
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NGC 3256 is the most luminous LIRG in the local uni¬ 
verse {z < 0.01), and it is a merging galaxy with two 
(N and S) nuclei. Presence of an AGN in the S nu¬ 
cleus has been controversial. We examined spectropho- 
tometric characteristics of the S nucleus at both near- 
mid-infrared (NIR-MIR) and X-ray using archival and 
published data, and found several pieces of evidence to 
support a low-luminosity AGN obscured by Compton- 
thin matter. The following are our findings and their 
implications. 

We found in IRAC flux ratio maps that the S nucleus 
shows distinct photometric properties at 3.6-8.0 /rm, 
most notably in its excess of 4.5 ^m flux. In contrast, the 
N nucleus is similar in colors to the star-forming regions 
within the host galaxy. We applied the IRAC color-color 
diagrams for AGN diagnostics to the nuclear photome¬ 
tries for each nucleus, and found that the N and S nuclei 
show starburst- and AGN (power-law)-like SEDs, respec¬ 
tively. This difference originates from the 4.5 /rm excess 
at the S nucleus. 

Using high-resolution HST NICMOS images, we ex¬ 
tracted a compact source at the S nucleus by subtracting 
the stellar component at 2.2 /rm. The S nucleus consists 
of an unresolved (at 0"26 FWHM) and a resolved (0"75 
FWHM or 128 pc after subtracting the instrumental res¬ 
olution) component. The corresponding structures are 
not seen at 1.6 /im and shorter wavelengths. Because 
of its position and size, we identify its unresolved core 
with the compact MIR core (at 0'.'30 FWHM resolution 
at 8.74 /im). The flux of the 2.2 /tm core is consistent 
with our AGN spectral energy distribution (SED) model. 

We analyzed the IRS nuclear spectrophotometries to¬ 
gether with the IRAC nuclear photometries at 3.6- 
14.5 /im. The IRS spectrum of the S nucleus shows bluer 
colors at < 6 /im with respect to both the N nucleus 
and the IRS starburst template of Brandi et al. (2006) 
in a consistent way with the IRAC SED. We conducted 
SED modeling of the S nucleus to reproduce the high 
spatial-resolution (0'.'36 aperture) T-ReCS Wband spec¬ 
trophotometry (Di'az-Santos et al. 2010a) as well as the 
nuclear IRAC 4.5 /im and IRS data. Our AGN-starburst 
composite model (a heavily absorbed power-law AGN 
SED superposed on a mildly absorbed starburst-powered 
LIRG SED template) successfully reproduces both the 
deep silicate absorption at 9.7 /tm within the 0'.'36 aper¬ 
ture and the distinct PAH features with the 4-6 /im ex¬ 
cess at ^ 3" scale. We estimated Ay toward the AGN to 
be as large as 80 mag. All the MIR results point toward 
a heavily absorbed (but Compton-thin) type-2 AGN in 
the S nucleus. The AGN is most likely accompanied by 
circumnuclear star-forming regions. 

We obtained a deep Chandra X-ray spectrum with 
~ 2.7 times more exposure time than in the previous 
study, and performed spectral analysis of a point-like 
source at the S nucleus. We found that a dual-component 
(for primary and scattered ones) power-law model suc¬ 
cessfully fits the hard spectrum. The inferred column 
density that is associated with the heavily absorbed pri¬ 
mary component is IVh = 7 ^3® x 10^^ cm“^, which is in 
the range typically observed in Seyfert 2s and is consis¬ 
tent with our MIR measurement. The lightly absorbed 
component can be interpreted as emission scattered by 
ionized medium in the opening part of the obscuring ma¬ 


terial around the AGN. The successful dual-component 
model suggests a Compton-thin type-2 AGN in the S 
nucleus. We also examined three other models for the S 
nucleus, namely a Compton-thick AGN, a single X-ray 
emitting compact source or a collection of such sources in 
the S nucleus, and thermal plasma with kT < 1 keV. We 
found that all of them are either rejected or less likely. A 
Compton-thick AGN is rejected because of a limit on EW 
of a fluorescent Fe-K emission line at 6.4 keV (< 190 eV 
for a reflection-dominated model). Models with compact 
stellar source(s) are less likely because of the observed 
high luminosity, hard spectrum, and/or small region size 
well within the Chandra resolution (0'.'49 or 83 pc), al¬ 
though we cannot reject possibilities of models with mul¬ 
tiple compact stellar sources. A model with the thermal 
plasma seems unusual, although it is statistically possi¬ 
ble, because the low-fcT plasma needs to be confined and 
hidden behind absorbing matter. 

Our model with AGN is quantitatively consistent with 
most of the earlier studies that found no AGN in the S 
nucleus if we consider different aperture sizes and the 
star-forming circumnuclear region that dominates the 
MIR flux at ~ 3" scale. In our best MIR SED model 
(the AGN-starburst composite model), the AGN con¬ 
tributes to ~ 2% and ~ 0.2% of the whole-aperture MIR 
luminosities at 6 and 24 /tm, respectively. In our best 
X-ray model (a dual absorber model with Seyfert-like 
power-law continua), the absorption corrected X-ray lu¬ 
minosity of the S nucleus is about 20% of the total X-ray 
luminosity integrated over NGG 3256. In particular, we 
showed that the X-ray to MIR luminosity ratio of the 
AGN is consistent with that of typical Seyfert 2s if we 
exclude contribution of MIR flux from the circumnuclear 
star formation. 
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Figure 1. HST NICMOS F160W and F222M images (upper row) and Spitzer IRAC four channel images (middle and lower rows) of 
NGC 3256 around the N and S nuclei. North is up and east is to the left. Both NICMOS and IRAC images are shown in a square root 
scale as indicated in scale bars (in units of MJy/steradian) for each panel. The N nucleus is saturated in the two NICMOS images, and 
are masked. Positions of the two nuclei in radio continuum (Neff et al. 2003) are marked with crosses on the NICMOS images. Contours 
at logarithmic steps are overlaid on the IRAC 3.6 pm image. In the lower right panel (IRAC 8.0 pm), two white circles indicate locations 
of artifacts (known as the “bandwidth effect”; IRAC instrument handbook) caused by the saturated N nucleus. Also in this panel, the 
saturated pixels around the N nucleus are fixed as described in the main text. 
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Figure 2. Flux ratio maps of NICMOS F222M/F160W (upper left), IRAC 4.5 /im/3.6 fim (upper right), 5.8 /rm/4.5 /rm (lower left), 
and 8.0 /im/5.8 /rm (lower right) around the N and S nuclei. All images are shown in the same sky area as in Figure 1. Contours of the 
IRAC 3.6 fim image shown in Figure 1 are overlaid on all flux ratio maps for positional references. Flux ratio scales are indicated in scale 
bars for each panel. Outer fainter regions with signal-to-noise ratio being less than 3.0 per pixel are masked in the F222M/F160W image. 
Two white circles in the 8.0 /j.m/5.8 /rm map (lower right) indicate areas that are affected by the artifacts in the 8.0 fim image as shown 
in Figure 1. 
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Figure 3. The IRS spectra and IRAC photometries of the N (black; top) and S (red; middle) nuclei, and the NICMOS photometry of the 
S nucleus (red). All photometries for the S nucleus are shown with an offset of —0.5 for clarity of the figure. IRAC fluxes of each nucleus 
are shown both for nuclear fluxes (diamonds connected by dot-dashed lines) and adjusted fluxes (to match the IRS spectra; diamonds 
connected by dashed lines. See the main text for the details of the adjustment. We also show the NICMOS F222M nuclear photometry of 
an unresolved compact source (lower) and sum of the unresolved and resolved sources (upper) at the S nucleus with red circles. The error 
bars are for dilcr. Each blue solid line shows the starburst IRS spectrum of Brandi et al. (2006) that is redshifted to NGC 3256 and scaled 
at 6.0-8.0 /im for each nucleus. The LIRG SED template of NGC 6090 arbitrary scaled for clarify of the figure is also shown with green 
solid line (bottom). 
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Figure 4. Zoom-up NICMOS images of the S nucleus and results of our analysis, (a)-(c) NICMOS FllOW, F160W, and F222M images, 
(d) Synthetic F222M image of the stellar component. See the main text for the details of the synthesis, (e) A flux ratio image of 
F222M/F160W. (f) The F222M nuclear component, made by subtracting the synthetic F222M stellar component image from the observed 
F222M image. The NICMOS resolution (0"26 FWHM) is indicated by a white circle, (g)-(h) The residual images of the F222M nuclear 
component after subtracting the fitted Gaussian models. The panel (g) shows a case of single-Gaussian model, and (h) shows a case of 
double-Gaussian model. The fitted positions of Gaussian components are marked with white pluses in panels (g) and (h). All panels are 
shown in linear scale, and flux (in units of MJy/steradian) and flux ratio (for only panel f) scale are indicated in scale bars for each panel. 
The zero flux levels in the residual images are also indicated by black contours. 
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Figure 5. (a) The spectrum of the flux ratio between the N and S nuclei. The IRS data are in black and the IRAC data are in blue, (b) 
The spectra of the flux ratio between the individual nuclei and the IRS starburst template of Brandi et al. (2006). The latter is redshifted 
and scaled at 6.0-8.0 /rm for each nucleus. The N nucleus is shown in black and the S nucleus is in red. Error bars are ±lcr. 
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Figure 6. The IRAC color-color diagrams of the N (open square) and S (open circle) nuclei with boundaries for AGN selection. Error 
bars are dilcr. Also shown in filled symbols are rest frame colors/flux ratios of ULIRG IRAS 22491-1808 from Polletta et ah (2007) (square), 
starburst M 82 from Polletta et ah (2007) (star), a template of normal star-forming spiral galaxy from Dale Sz Helou (2002) (triangle), and 
a representative elliptical galaxy from Polletta et ah (2007) (circle). Connected diamonds are for pure power-law SEDs, and each triangle 
corresponds to a power-law index of a from —0.5 (lower left) to —3.0 (upper right) in steps of 0.5. Connected pluses represent colors/fiux 
ratios for pure blackbody SEDs, and each plus corresponds to a blackbody temperature of 600 K (upper right), 1000 K, 3000 K, 6000 K, 
and 10000 K (lower left). The boundaries for AGN selection are from Donley et ah (2012). 
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Figure 7. Comparison of extinction curves used in our SED modeling. We show the extinction curve of Chiar Sz Tielens (2006) for the 
Galactic center (black solid line), the original PAHFIT (Smith et ah 2007b) extinction curve (blue solid line), the original extinction curve 
of Fritz et ah (2011) (broken magenta line), and the modified PAHFIT extinction curve (red solid line) to reproduce the curve of Fritz et 
ah (2011). See the main text for the details of the modified PAHFIT extinction curve. 
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Figure 8. SED modeling results for the S nucleus. The absorbed single-starburst model (Single-SB), the AGN-starburst composite model 
(AGN-hSB), and the starburst-starburst composite model (SB-I-SB) are shown in the top, middle, and bottom rows, respectively. Results 
with the extinction curve of Chiar Sz Tielens (2006) and the modified PAHFIT extinction curve are shown in the left and right columns, 
respectively. The adjusted IRAC photometries as shown in Figure 3 (connected red diamonds with ±lcr error bars), the IRS (with dilcr 
error bars; black) and T-ReGS spectrophotometries (black) are shown in all panels. For the absorbed single-starburst model, we overlay 
the best-fit models with magenta solid lines. For both the AGN-starburst and starburst-starburst composite models, we overlay the best 
fit models for the T-ReGS (blue), residual 4.5 /im-IRS spectrophotometries (observed—best-fit T-ReGS model) (cyan), and sum of the two 
components (magenta). The sigma ((observation—model)/observation noise) spectra of the IRS spectrophotometries (black) and synthetic 
IRAC fluxes from the best-fit models (red diamonds) are shown for each panel. 
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Figure 9. The best SED modeling result for the S nucleus with the AGN-starburst composite model and the modified PAHFIT extinction 
curve at 2.0-10 pm. The same colors and symbols of Figure 8 are used. The NICMOS photometry at 2.2 pm of the unresolved core is also 
shown with a dilcr error bar (red circle). 













Energy (keV) 

Figure 10. Observed Chandra spectra of the S nucleus for the first (open circles) and second (filled circles) observations. Data are binned 
for presentation purpose. (Upper panel) Observed spectra and the best-fit model (Model D, solid line). Only the model for the second 
observation is shown for clarity. (Lower panel) Data/Model ratio. 

The error bars are for dilcr. 
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Figure 11. Unfolded Chandra spectrum of the S nucleus. The responses of detector and telescope are unfolded. The spectra from the two 
observations are combined for presentation purpose. The error bars are for ±lcr. The best-fit model (Model D, solid line), lightly (dashed 
line), and heavily (dotted line) absorbed components are shown for comparison. 





26 


Ohyama, Terashima, & Sakamoto 


Table 1 

IRAC photometry results. 


Channel 

4'lux (mjy)®' 


N nucleus 

S nucleus 

Nuclear photometry*^ 

3.6 /im 

17.0 ±0.2 

9.4 ±0.1 

4.5 p-m 

17.9 ±0.2 

24.5 ±0.2 

5.8 p-m 

63.5 ±0.9 

41.0 ±0.6 

8.0 fim. 

203 ± 22'= 

51.0 ±2.7 

Adjusted for the IRS spectra*^ 

3.6 /im 

34.7 ±0.4 

25.2 ±0.2 

4.5 p-m 

36.6 ±0.4 

41.2 ±0.3 

5.8 fim 

130 ± 1.9 

100 ± 1.1 

8.0 fim 

415 ± 44" 

240 ± 20 


^Uncertainties are for dilcr. 

^Photometry with 2'!% diameter aperture with aperture correction. 

^The photometry is made on the image corrected for the saturation at the N nucleus. See the main text for the details. 

"^IRAC fluxes adjusted for the same level of circumnuclear flux contamination for the IRS resolution. See the main text for the details of 
the adjustment. 


Table 2 

Results of the double-Gaussian model fitting to the S nucleus on the stellar-component subtracted F222M image. 



compact source extended source 

Position 

(J2000) 

RA 

10;27;51.24=- 

DEC 

-43:54:19.2'' 

source size 

(FWHM in arcsec) 
unresolved^ 0'.'75 ± 0701 " 

Flux 

(mjy)'^ 

0.65 ±0.02 5.33 ±0.16 


^Uncertainties (0'^2) are dominated by global astrometry solution. See the main text for the details. 
^Assumed to be unresolved at 0^^26 FWHM. 

^After correcting for the instrumental resolution of 0^.^26 FWHM. 

^Uncertainties are for dilcr. 


Table 3 

MIR spectrophotometry model fitting results. 


Extinction curve 

Data to fit 

models^ 

Single-SB" 

AGN-hSB 

(y/36 aperture ~ 3'' aperture 

SB+SB 

0^/36 aperture ~ S" aperture 

Chiar & Tielens (2006)" 

4.5 /rm-l-lRS-l-T-ReCS 

'^9.7 

1.04 ±0.01 

9.4 ±0.4 

0.68 ±0.01 

6.8 ±0.3 

0.98 ±0.02 



reduced 

4.09 

1.4 

1.4 

1.2 

4.4 


IRS+T-ReCS 

reduced x^ 

0.90 

1.4 

0.77 

1.2 

0.87 

Original PAHFIT 

4.5 /rm-l-IRS+T-ReCS 

^9.7 

0.00 

8.0 ±0.4 

0.0 

6.0 ±0.3 

0.0 



reduced x^ 

5.80 

3.6 

5.2 

2.7 

6.0 


IRS+T-ReCS 

reduced x^ 

1.09 

3.6 

1.08 

2.7 

1.11 

Modified PAHFIT'^ 

4.5 /rm-l-IRS+T-ReCS 

^9.7 

1.21 ± 0.05 

12.7 ±0.5 

0.48 ± 0.04 

8.3 ±0.4 

1.10 ±0.02 



reduced x^ 

3.58 

1.2 

0.79 

1.2 

3.5 


IRS+T-ReCS 

reduced x^ 

0.90 

1.2 

0.74 

1.2 

0.86 


^Models: Single-SB = absorbed single-starburst model; AGN-hSB = AGN-starburst composite model; SB-l-SB = starburst-starburst 
composite model. 

^T-ReGS data are not used for this model. 

^Toward Galactic center. 

^Modified to reproduce the extinction curve of Fritz et al. (2011). See the main text for the details. 
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Table 4 

Chandra Observation Log. 


Observation ID^ 

Date‘s 

Exposure^ 

(ksec) 

Counts'' 

835 

2000 Jan 5 

16.2 

39 

3569 

2003 May 23 

27.2 

62 


^Observation identification number. 

'^Observation start date. 

‘^Exposure time after data screening. 

‘'Number of counts after background subtraction in 0.5-7 keV. 


Table 5 

Results of X-ray spectral fits with power law model. 


ModeU 

Ah.i 

Ah .2 / 

Photon index 

Norm. 

Constant 

C/dof 

EW 

F2-IO keV L 

2-10 keV 


(10^^ cm 

(10^^ cm 





(eV) 



(1) 

(2) 

(3) (4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

A 

0.054 (<0.97) 


U.4Z_0.39 

0.067 

0.94t°;« 

76.8/90 

<150 

11.5 

1.7 

B 

3.6 


1.8'" 

2.15 

0.92 

94.5/91'’ 




C 

0 (< 0.27) 


1.8'" 

46.4 

0.93l°'3? 

78.1/91 

<190 

9.26 

1.4 

D 

U.o4_q 25 

7+iy Q qho+u.u; 5 y 

/_3 u.y4Y)_Q ^49 

1.8'" 

2.93 

0.95l°;^° 

78.6/89 

< 550 

6.33 

1.5 

Note. 

— (1) Spectral models. The Galactic absorption (Ah 

= 9.14 X 

: 10^° cm’ 

^) is applied to all the models. 

(2), (3) 


Intrinsic absorption column density. (4) The fraction of continuum absorbed by a large column density (Ah, 2). (5) Photon 
index of power law component. (6) Normalization of power law in units of 10~® photons cm“^ at 1 keV. (7) Constant factor 
relative to the first observation. (8) C statistic and degrees of freedom. (9) Upper limit on the equivalent width of fluorescent 
Fe-Ka line in units of eV. (10) Observed flux in the 2-10 keV band for the first observation in units of 10“^^ erg s“^ cm“^. 
(11) Luminosity corrected for absorption in the 2-10 keV band for the first observation in units of 10^° erg s~^. 

Fixed parameter. 

Errors, limit on Fe-K line equivalent width, flux, and luminosity are not shown for this model since the value of C statistic 
is significantly worse than those for other models. 


Table 6 

Results of X-ray spectral fits with thermal model. 


ModeU 

(1) 

Ah 

(1022 cm-2) 
(2) 

Ah 

(1022 cm-2) 

(3) 

kT 

(4) 

Norm. 

(5) 

kT 

(6) 

Norm. 

(7) 

Constant 

(8) 

c/dof 

(9) 

F2_10 keV 

(10) 

L 2 -IO keV 

(11) 

E 

F 

K 1 +3.3 
^•-*--1.9 

4 3+4.0 

1 3+1.0 
-‘-•'^-0.5 

64“ (>9.4) 
64“ (>9.7) 

6.71 

7.44 

64“ (>1.2) 
1.0*' 

0.428 

0.798 

u.yo_Q 26 

0 95+0-40 

u.yo_Q 28 

77.6/88 

77.3/88 

6.98 

7.06 

1.4 

1.4 


Note. — (1) Spectral models. The Galactic absorption (Ah = 9.14x 10'^° cm“^) is applied to all the models. (2) Intrinsic 
absorption column density for heavily absorbed APEC component. (3) Intrinsic absorption column density common to 
both APEC components. (4) Temperature of heavily absorbed APEC component. (5) Normalization of heavily absorbed 
APEC component in units of 10“^®/(47r(L>A(l+ 2))^//UeUndU, where Da is the angular size distance to the source (cm), 
rie is the electron density (cm~®), and nn is the hydrogen density (cm~®). (6) Temperature of lightly absorbed APEC 
component. (7) Normalization of lightly absorbed APEC component in units of 10“^®/(47r(HA(l + 2))^// UenudV. (8) 
Constant factor. (9) C statistic and degrees of freedom. (10) Observed flux in the 2-10 keV band for the first observation 
in units of 10”^^ erg s~^ cm“^. (11) Luminosity corrected for absorption in the 2-10 keV band for the first observation in 
units of 10"^° erg s“^ 

a Pegged at the highest temperature allowed in the fit. 

Fixed parameter. 







